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1.  INTRODUCTION 


1.1  Object  and  Scope 

Wave  propagation  problems  in  linearly  elastic,  homogeneous, 
anisotropic  media  have  attracted  increased  interest  in  recent  years. 

Aside  from  the  inherent  interest  in  the  subject  as  part  of  theoretical 
mechanics,  results  of  investigations  in  this  area  have  sufficient  practical 
engineering  applications  to  motivate  further  explorations  of  the  math¬ 
ematical  and  computational  difficulties  involved.  In  the  past,  the  effects 
of  anisotropy  were  often  neglected  as  being  of  secondary  importance  in 
the  dynamics  of  elastic  media.  However,  when  detailed  computations  are 
carried  out  for  wave  propagation  in  isotropic  media  using  large  digital 
computers  and  involving  massive  arithmetic  effort,  it  would  seem  that 
more  attention  should  be  given  to  what  are  possibly  significant  effects  of 
anisotropy. 

Several  theoretical  investigations  have  been  stimulated  by  the 
development  of  ultrasonic  techniques  for  the  measurement  of  the  dynamic 
elastic  constants  in  pure  crystals,  where  anisotropy  plays  a  fundamental 
role.  These  techniques  were  based  on  the  availability  of  the  necessary 
electronic  equipment  to  excite  and  detect  mechanical  vibrations  of  high 
frequency  which  enabled  experimental  researchers  to  observe  the 
characteristics  of  propagation  of  disturbances  with  wave  lengths  very 
small  in  comparison  with  the  dimensions  of  the  source  and  specimen.  In 
theoretical  seismology,  increased  attention  has  been  given  to  studies 
concerning  the  effects  of  a  possible  continental  anisotropy  in  the 
propagation  of  seismic  waves,  in  particular  the  Rayleigh  surface  waves. 
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Seismic  waves  are  also  of  great  interest  in  earthquake  engineering. 
Specifically,  the  effect  of  anisotropy  on  strong  ground-motion  records 
requires  further  exploration.  The  propagation  of  waves  in  piezoelectric 
crystals,  the  dynamic  behavior  of  composite  materials,  and  a  variety  of 
problems  in  solid  state  physics  and  geophysics  fall  into  the  range  of 
applications  where  an  increased  knowledge  of  dynamics  of  anisotropic 
materials  is  desirable. 

The  object  of  the  present  study  is  to  provide  analytical  and 
computational  techniques  to  obtain  solutions  for  the  transmission  of 
waves  from  a  point  force  in  a  linearly  elastic,  homogeneous,  anisotropic 
medium.  Two  types  of  problems  have  been  undertaken.  The  first  concerns 
the  two-dimensional  problem  of  a  line  force  applied  to  the  surface  or 
in  the  interior  of  an  orthotropic  material,  i.e.  a  material  which  has  at 
each  point  three  orthogonal  planes  of  elastic  symmetry.  The  boundary 
plane  is  taken  parallel  to  one  of  the  symmetry  planes  and  the  applied 
load  is  a  step  function  in  time.  The  stresses  are  determined  directly 
as  functions  of  a  complex  variable  while  the  displacements  require  the 
evaluation  of  a  simple  quadrature.  The  second  type  of  problems  considered 
in  this  study  concerns  a  step  force  applied  to  the  surface  of  a 
transversely  isotropic  material.  Transverse  isotropy  is  a  special  case 
of  anisotropy  in  which  there  is  at  every  point  of  the  medium  one  axis 
of  elastic  symmetry  while  all  directions  perpendicular  to  that  axis  are 
equivalent.  The  dynamic  response  given  by  the  three  displacements  in  a 
cylindrical  coordinate  system  requires  the  calculation  of  a  simple 
quadrature  in  the  complex  plane. 


1.2  Previous  Related  Studies 


Investigations  concerning  anisotropic  wave  propagation  can  be  traced 
back  to  the  research  work  on  crystal  optics  at  the  end  of  the  last  century 
which  laid  down  the  foundation  for  a  qualitative  understanding  of  the 
effects  of  anisotropy  on  elastic  waves  propagating  in  a  solid  medium. 
Rigorous  theoretical  developments  awaited  the  availability  of  ultrasonic 
techniques  which  occured  several  decades  later  and  spurred  a  renewed 
interest  in  the  subject. 

The  basic  principles  of  the  theory  of  elastic  wave  propagation  in 

* 

anisotropic  media  can  be  found  in  the  monographs  by  Hearmon  [15]  , 

Fedorov  [13],  Musgrave  [25],  Auld  [3],  and  Chadwick  and  Smith  [8]. 

Recent  advances  of  elastic  waves  at  crystalline  interfaces  were  discussed 
by  Musgrave  [26]  who  also  included  a  list  of  original  contributions  to 
various  specialized  problems.  Kraut  [19]  examined  the  problem  of  a 
vertical  line  force  on  the  surface  of  a  transversely  isotropic  material, 

De  [10]  discussed  the  Rayleigh  wave  contribution  to  the  displacements 
produced  by  a  suddenly  applied  force  at  the  surface,  while  Payton  [34,  35] 
gave  the  solutions  for  suddenly  applied  point  force  and  for  time  dependent 
line  load  in  an  unbounded  transversely  isotropic  medium.  In  all  these 
solutions  the  method  of  analysis  is  based  on  the  classical  multiple 
transform  techniques,  with  the  attendant  computational  difficulties. 

A  different  approach  was  followed  by  Sveklo  [44]  who  examined  the 
two-dimensional  problem  of  an  instantaneous  pulse  in  an  anisotropic 
medium  by  employing  the  Smirnov-Sobolev  method  of  complex  solutions,  which 

ifr 

Numbers  in  brackets  refer  to  entries  in  the  List  of  References. 
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was  initially  applied  for  Isotropic  solids  [39,  40].  Sveklo  used  the 
same  approach  to  solve  the  problem  of  a  force  applied  on  the  surface  of 
a  half-space  (Lamb's  problem)  [45].  The  anisotropy  of  the  material  was 
defined  by  three  elastic  constants.  The  same  type  of  material  was  used 
by  the  same  author  to  examine  dynamic  problems  with  mixed  boundary 
conditions  [46].  Reflected  and  refracted  waves  at  the  interface  of  two 
anisotropic  media  with  three  elastic  constants  were  examined  by  Osipov 
[27]  who  also  in  a  number  of  other  publications  covered  topics  such  as 
the  behavior  of  the  wave  propagation  velocities  in  materials  with  four 
elastic  constants  [28],  the  application  of  the  Smirnov-Sobolev  method  to 
formulate  solutions  for  anisotropic  materials  [29,30],  the  examination 
of  the  Rayleigh  type  waves  [31],  and  some  cases  of  mapping  of  the  wave 
fields  on  the  complex  plane  [32].  A  number  of  errors  have  been  noted  in 
Osipov's  work;  some  of  them  were  recognized  and  corrected  by  the  author 
himself  in  later  publications.  Lamb's  problem  in  anisotropic  media  was 
also  examined  by  Budaev  [5,  6,  7]  in  terms  of  certain  dimensionless 
quantities  derived  from  the  elastic  constants. 

1.3  Method  of  Solution 

Two-dimensional  wave  propagation  problems  in  homogeneous  elastic 
anisotropic  bodies  require  the  solution  of  a  system  of  two  partial 
differential  equations.  The  Smirnov-Sobolev  method  of  self-similar 
potentials  utilizes  some  important  results  of  complex  analysis  to  derive 
the  solutions  as  analytic  functions  of  a  complex  variable.  The  complex 
variable  is  determined  from  an  equation  which  describes  the  general 
solution  of  the  equations  of  motion  in  terms  of  the  characteristic 
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surfaces. 

A  given  set  of  boundary  conditions  which  are  homogeneous  functions 
of  the  x,  y  and  t  variables  provides  the  necessary  equations  to  determine 
the  analytic  functions  or  potentials  mentioned  previously.  The  real  part 
of  the  complex  solutions  constitutes  the  answer  to  a  problem  because  the 
real  part  of  the  complex  stresses  and/or  displacements  represents  the 
prescribed  boundary  conditions. 

The  anisotropy  of  the  medium  in  the  two-dimensional  problems  is 
defined  by  four  elastic  constants.  The  applied  load  is  given  as  a  line 
force  parallel  to  the  elastic  direction  so  that  the  conditions  of  plane 
strain  are  satisfied.  When  the  medium  is  unbounded,  two  types  of  waves 
propagate  from  the  source  of  disturbance.  The  first  type  is  called 
quasi-longitudinal  and  the  second  quasi -transverse  in  analogy  with  the 
isotropic  case.  The  existence  of  a  plane  boundary  introduces  reflected, 
surface  and  head  waves. 

The  solution  of  the  three-dimensional  equations  of  motion  for  a 
transversely  isotropic  medium  is  reduced,  by  the  technique  of  rotational 
superposition  described  In  Chapter  3,  to  the  solution  of  equations  of 
motion  for  a  plane  strain  problem  and  an  antiplane  problem.  Both  of 
these  cases  are  treated  in  Chapter  2.  Thus,  the  simplicity  of  the  method 
of  self-similar  potentials  is  preserved  in  developing  three-dimensional 
solutions. 

The  study  is  organized  as  follows.  Chapter  2  briefly  explains  the 
Smirnov-Sobolev  method  as  it  is  applicable  to  anisotropic  materials  of  a 
certain  kind  and  gives  the  solutions  to  a  variety  of  two-dimensional 
problems.  Chapter  3  explains  the  technique  of  rotational  superposition 


and  presents  the  solution  of  certain  axisymmetric  and  non- axl symmetric 
problems  in  the  form  of  the  displacements.  The  often  intricate  mapping 
of  the  wave  fields  in  the  complex  domain  defined  by  a  four-sheeted 
Riemann  surface  is  the  subject  of  Chapter  4.  In  Chapter  5  numerical 
results  are  discussed  together  with  some  of  the  techniques  required 
to  overcome  inherent  computational  difficulties.  A  summary  of  the  study, 
conclusions  and  recommendations  for  further  research  are  given  as 
Chapter  6. 

1.4  Notation 


The  symbols  used  in  this  study  are  defined  in  the  text  where  they 


first  appear. 

For  convenient  reference,  those  most  frequently  used  are 

listed  below. 

d)  by  C j  d  y  6 

elastic  parameters  of  orthotropic 
and  transversely  isotropic  materials 

Ay  By  Cy  Dy  E 

the  five  material  cases  according  to 
the  configuration  of  the  quasi- 
transverse  wave  front 

V  Be’  ce’  De 

curves  on  the  complex  domain  which 
separate  the  regions  of  mapping  the 
quasi-longitudinal  and  the  quasi- 
transverse  wave  field 

A01 ’  a02 

points  in  the  wave  fields  which  map 
on  the  branch  points  Oq-j  and 

cij  J  =  2»  •••6) 

CK  (k  =  1,  2,  3) 

EX 


elastic  constants  of  a  general 
anisotropic  material 

contours  of  integration  in  the  three- 
dimensional  problems 


normal  strain  in  the  X-direction 
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weighting  functions 

magnitudes  of  horizontal  and  vertical 
forces 

the  plane  strain  problem  of  a 
horizontal  step  force  in  the  full- 
space 

function  defined  in  Eq.  (3.25) 

the  plane  strain  problem  of  a  vertical 
step  force  on  the  surface  of  a  half¬ 
space 

the  plane  strain  problem  of  a 
horizontal  step  force  in  the  interior 
of  a  half-space 

quantities  involving  the  elastic 
parameters,  defined  in  Eq.  (2.10) 

functions  defined  in  Eq.  (2.24) 

unit  vector  normal  to  the  tangent 
at  the  wave  front 

quantities  involving  the  elastic 
parameters,  defined  in  Eq.  (2.10) 
and  (C.10) 

function  defined  in  Eq.  (2.20) 

the  radial  coordinate  in  the 
cylindrical  system 

Rayleigh  function  in  anisotropic 
media  defined  in  Eq.  (2.28) 

functions  defined  in  Eqs.  (2.47) 
and  (2.48) 

the  four  sheets  of  the  complete 
Riemann  surface 

function  defined  in  Eq.  (2.24) 


the  time  variable 
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arrival  times  of  the  wave  fronts 


the  times  of  the  beginning  of  the 
mapping  of  the  two  wave  fields  in 
the  complex  domain 

displacement  components  of  the  plane 
strain  and  the  antiplane  problems 
(an  asterisk  indicates  the  complex 
quantities) 

complex  displacements  at  the  free 
surface  (y  =  0) 

displacement  components  in  the  three- 
dimensional  problems  (an  asterisk 
indicates  the  complex  quantities) 

displacement  components  of  the  plane 
strain  and  the  anti  plane  problems 
as  used  to  develop  the  three- 
dimensional  solutions 

angular  quantity  (=  u  -  w0) 

the  cartesian  coordinate  system 

coordinates  of  a  point  in  the  (x,  y) 
plane 

a  cartesian  coordinate  system  rotated 
by  wq  about  the  y  =  Y  axis 

the  y-coordinate  of  the  point  where 
the  force  is  applied  in  the  interior 
of  the  half-space 

the  angular  quantities  shown  in  Fig.  2 

values  of  the  quantity  bk  corresponding 
to  various  values  of  e  (see  Table  1) 

~  XR^R 
shear  strains 

a  value  of  y  such  that  the  point 
(xr>  yR)  is  located  inside  the  sector 
formed  from  the  cusps  and  the  origin 
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r 

o 

6k  (k  =  1,  2,  3) 

&'K  (K  =  1,  2,  3) 

A^,  A^  (k=  1,  2 ,  3) 


n 

nl  >  n2>  nj »  n2 

*  1*  2* 
n»  n  ,  n 


0 

1*  2*  *  **  '1*  ~2* 

0 1 ,  e  ,  e ,  e  ,  e  ,  e  ,  0 


On,  o 


Or 


D’  wor  u02 


K 


Ak(0)  (<  =  1,  2,  3) 

5 


any  value  of  y  that  is  not  rc 

function  defined  by  Eqs.  (2.16) 
and  (2.31) 

=  96  / 3 0 

K  1C 

functions  equivalent  to  6K,  6^, 
as  used  for  the  three-dimensional 
problems 

real  and  imaginary  parts  of  0^ 


normal  strains 


real  and  imaginary  parts  of  *K(0K) 

dimensionless  coordinate  (=  y/t) 

values  of  n  at  the  crossings  of 
the  wave  fronts  with  the  n-axis 

points  on  the  n-axis  which  correspond 
to  various  values  of  0  according  to 
conditions  given  in  Chapter  4 

complex  parameter 

real  values  of  0  corresponding  to 
cusps  and  nodes  (double  points)  of 
the  quasi-transverse  wave  front 
(see  Table  1 ) 

the  branch  points  on  the  Riemann 
surface 

index  referring  to  the  quasi - 
longitudinal  wave  (<=1),  the  quasi- 
transverse  wave  (<=2)  or  the  shear 
wave  in  the  antiplane  problem  (<=3) 

function  defined  in  Eqs.  (2.19)  and 
(2.32) 


S2,  C],  q 


dimensionless  coordinate  (=  x/t) 

values  of  j;  at  the  crossings  of  the 
wave  fronts  with  the  c-axis 


the  ^-coordinate  of  a  node  located  on 
the  5-axis 

function  defined  in  Eq.  (2.27) 


mass  density  of  the  material 

normal  stresses  in  the  two- 
dimensional  problems  (an  asterisk) 
indicates  the  complex  quantities) 

complex  normal  stress  at  the  free 
surface  (y=0) 

normal  stresses  in  the  cylindrical 
coordinate  system 

summation  symbol 


normal  stresses  for  plane  strain 
problems  as  used  to  develop  the 
three-dimensional  solutions 

complex  self-similar  fictitious 
normal  stress 

shear  stresses  of  the  plane  strain 
and  the  anti  plane  problem  (an 
asterisk  indicates  the  complex 
quantities) 

shear  stresses  in  the  cylindrical 
coordinate  system 

shear  stresses  of  the  plane  strain 
and  anti  plane  problems  as  used  to 
develop  the  three-dimensional  solutions 

complex  self-similar  fictitious 
shear  stresses 

variable  used  in  Chapter  3  (=  rcosW) 

functions  defined  in  Eqs.  (4.9) 
and  (4.12) 


functions  defined  in  Eqs.  (4.4)  and 
(4.11) 
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t  Jim 


wmamm 


m* 


\(0)»  «Kj(e) 


n 


the  angular  coordinate 
cylindrical  system 

self-similar  potential 


in  the 

functions 


2.  TWO-DIMENSIONAL  WAVE  PROPAGATION  PROBLEMS 
IN  ANISOTROPIC  MATERIALS 

2.1  Generalized  Hooke's  Law 

In  continuous  media  the  stresses  and  strains  at  a  point  are  completely 
expressed  by  the  corresponding  tensors.  Each  component  of  the  stress 
tensor  can  be  written  as  a  function  of  the  components  of  the  strain  tensor. 
Under  the  assumptions  that  the  deformation  takes  place  at  a  fixed 
temperature,  and  that  the  initial  unstrained  state  is  also  unstressed,  the 
linear  parts  of  these  functional  relations  between  stresses  and  strains 
comprise  a  system  of  equations  known  as  the  generalized  Hooke's  law. 

In  a  cartesian  coordinate  system,  where  the  engineering  notation  for 
the  stresses  and  strains  is  used,  the  generalized  Hooke's  law  can  be 
written  in  the  form  [43] 

°x  =  C1 1 Gx+Cl 2Gy+cl 3Gz+c14Yyz+cl 5Yzx+c16Yxy 
“y  ’  c21Gx+c22Gy+c23Gz+c24Yyz+c25Yzx+c26Yxy 
°z  *  c31Gx+c32Gy+c33Gz+c34Yyz+c35Yzx+c36Yxy 

Tyz  ■  c41Gx+c42Gy+c43Gz+c44Yyz+c45Yzx+c46Yxy  (2.1 ) 

Tzx  * c51Gx+c52Gy+c53Gz+c54Yyz+c55Yzx+c56Yxy 
Txy  =  c61Gx+c62Gy+c63Gz+c64Yyz+c65Yzx+c66Yxy 

The  coefficients  c^.  (i,  j  =  1,  2,  ...6)  are  the  elastic  constants  or 
moduli  of  the  material.  In  homogeneous  materials  they  are  independent 
of  the  system  coordinates.  Moreover,  the  existence  of  a  quadratic  strain 
energy  function  imposes  the  symmetry  conditions 
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cij  =  (i,  j  =  1,  2,  ...6)  (2.2) 

In  the  most  general  case  of  anisotropic  materials,  the  number  of 

independent  coefficients  required  to  determine  the  elastic  behavior  is 
★ 

twenty-one.  The  existence  of  axes  or  planes  of  elastic  symmetry  together 
with  an  appropriate  orientation  of  the  coordinate  system  such  that  some 
or  all  of  the  coordinate  axes  are  parallel  to  the  symmetry  axes  or  planes 
causes  a  reduction  of  the  number  of  constants  required  to  describe  the 
material.  The  highest  degree  of  elastic  symmetry  is  exhibited  by 
isotropic  media,  which  require  only  two  elastic  constants.  A  discussion 
of  the  method  used  to  obtain  the  necessary  constraint  equations  among 
the  constants  c..,  which  must  hold  whenever  various  degrees  of  elastic 

I  J 

symmetry  are  present,  can  be  found  in  Love  [21,  §105]  and  Voigt  [49]. 

2.2  Plane  Strain  Problems 

The  solution  of  dynamic  problems  in  anisotropic  materials  with 
triclinic  symmetry,  i.e.  involving  all  twenty-one  elastic  constants, 
presents  formidable  computational  difficulties,  mostly  as  a  result  of  the 
large  number  of  parameters  involved.  However,  great  insight  into  the 
effects  of  anisotropy  can  be  gained  by  the  study  of  materials  which 
possess  higher  degrees  of  symmetry. 

2.2.1  Orthotropic  Materials 

Consider  an  elastic  anisotropic  medium  which  at  each  point  has  three 
mutually  perpendicular  planes  of  elastic  symmetry.  These  may  be  taken 
* 

This  number  is  based  on  the  energy  approach  of  Green.  See  Love  [21]  for 
a  discussion  of  the  controversy  between  the  followers  of  Green's  approach 
and  the  "atomic"  modelers  who  advocated  the  "rari-constant"  theory  using 
fewer  constants. 
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parallel  to  the  axes  of  an  orthogonal  coordinate  system.  In  this  system 
it  can  be  shown  that  the  following  coefficients  vanish 

C14=C15=C16=C24=C25=C26=C34=C35=C36=C45=C46=C56=0  ^2‘3^ 

The  remaining  nine  constants  characterize  an  important  class  of  anisotropic 
materials  most  widely  known  as  orthotropic.  When  the  conditions  of  plane 
strain  are  fulfilled,  i.e. 


u  =  y  =  y  =  0 
2  yz  zx 


(2.4) 


the  stress-strain  relations  (2.1)  are  written  in  the  simple  form 


°X  "  cll'x+c12ey 
°y  "  c12cx+c22cy 
xy  66  xy 
The  other  stresses  are 

T  =  T  =0 
yz  zx 

CTz  =  c13ex+c23ey 


(2.5) 


(2.6) 


Equations  (2.5)  involve  four  independent  elastic  constants  and  determine 
the  plane  strain  problem  of  orthotropic  materials  in  the  x-y  coordinates.* 
In  the  dynamic  problems  of  this  study,  it  is  convenient  to  use  a  new  set 
of  parameters,  indicated  by  the  letters  a,  b,  c  and  d  which  are  defined  as 


*  ,  ... 

The  elastic  constants  c-| o  and  c?3  are  not  in  any  way  involved  in  the 

solution  of  the  plane  problem  ana  are  required  only  for  determining  the 

stress  a  normal  to  the  plane  of  deformation(x,  y). 
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T  C11  ,  _  C22  _  c12+c66  ,  _  C66 
a  ‘  ~  b  "  ~  C - p - *  d  "  ~ 


(2.7) 


where  p  is  the  mass  density  of  the  material.  These  parameters  have  units 

O  O 

of  velocity  to  the  second  power  [L  /Tc]  and  are  subject  to  the  conditions 


a  >  0,  b  >  0,  d  >  0,  ab-(c-d)  >0 


(2.8) 


stemming  from  the  requirement  of  a  positive  definite  strain  energy 
function. 

It  has  been  observed  that  for  all  known  media  the  constraints 


a  >  d  and 


b  >  d 


(2.9) 


although  not  necessary  for  the  positive  definite  character  of  the  strain 

★ 

energy,  seem  to  be  satisfied. 

The  following  quantities  involving  the  elastic  parameters  will 
appear  frequently  throughout  this  study  and  require  an  early  definition 

2  2 

L  =  ab  +  d  -c^ 


K.j  =  ab  -  (c-dr 
K2  =  ab  -  (c+d)2 
N1  =  (a-d)  (b-d)  -  c‘ 
N2  =  b(a-d)  -  c2 


(2.10a) 


This  is  reminiscent  of  the  better  known  fact  that  values  of  Poisson's 
ratio  in  the  range  (0,  -1)  do  not  violate  the  condition  for  a  positive 
definite  strain  energy,  but  are  never  observed  in  isotropic  materials. 


N4  =  d2(b-d)2  -  K1  K2  (2.10b) 

M  =  (b+d)  N1-d(a-b)  (b-d)  =  b^+d  [(b-d)2  -  c2] 


It  will  be  seen  later  that  the  range  of  values  assumed  by  certain 
of  these  quantities  greatly  influences  the  procedure  for  computation  of 
the  wave  fields,  even  requiring  separation  into  different  cases  or  sub¬ 
classes  of  orthotropic  materials. 


I 

2.2.2  The  Equations  of  Motion 


The  consideration  of  the  equilibrium  of  an  infinitesimal  element  in 
the  two-dimensional  space  and  the  application  of  the  linearized  strain- 
displacement  formulas 
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The  displacement  vector  u  (ux> 
formula  [37] 


Uy)  can  be  decomposed  by  Helmholtz's 


u  =  grad  4  +  curl  ip  (2.13) 

into  a  sum  of  two  vectors.  The  first  vector  represents  the  irrotational 
or  longitudinal  part  given  from  a  scalar  potential  function  <p,  while  the 
second  is  formed  from  a  vector  potential  function  ip,  satisfying  the 
condition  div  ip  =  0,  and  represents  the  equivoluminal  or  tran verse  part 
of  u.  Equation  (2.13)  applied  to  two-dimensional  wave  propagation  problems 
is  equivalent  to  the  following  two  relations 


=  li  +  ii 

ux  3x  3y* 


.it  14 

y  3y  3x 


(2.14) 


where  ip  stands  for  the  z-component  of  the  vector  ip. 

Direct  application  of  these  last  expressions  to  the  equations  of 
motion  (2.12)  leads  to  a  system  of  equations  in  terms  of  the  potentials 
<p  and  ip,  i  .e. 


JL 

3x 

3 


3y 


(a-c)  h  +  d  it 
3X^  sy 

(b-c)  &  - 
3x  3y 


=  0 

(2.15) 


S2<P  1 

at2  J 


=  0 


It  should  be  noted  that  the  equations  are  coupled,  unlike  the  isotropic 
case  where  the  same  process  leads  to  two  equations  each  involving  only  one 
of  the  potentials.  The  physical  meaning  is  that,  in  general,  each  of  the 
two  solutions  occurring  in  anisotropic  media  includes  both  irrotational 
and  equivoluminal  components  in  the  displacement  vector. 


•  J  ' 


2.3  The  Method  of  Self-Similar  Potentials 


The  method  of  self-similar  potentials  or  functionally  invariable 
solutions  is  a  powerful  technique  developed  by  the  Soviet  mathematicians 
V.  I.  Smirnov  and  S.  L.  Sobolev  [39,  40]  to  solve  the  equations  of  motion 
in  problems  of  two-dimensional  elasticity  with  initial  or  boundary 
conditions  given  as  homogeneous  functions  of  the  spatial  and  time  variables. 
An  extended  description  of  the  method  with  applications  to  certain  contact 
problems  is  given  by  Thompson  and  Robinson  [47].  The  application  of  the 
same  method  in  solving  dynamic  problems  in  a  three-dimensional  space  is 
given  by  Johnson  and  Robinson  [18]  and  Farewell  and  Robinson  [12],  while 
Seyyedian  and  Robinson  [38]  solved  the  problem  of  the  buried  dislocation 
pulse  in  a  layered  half  space  by  a  combination  of  self-similar  solutions. 

In  all  these  cases  the  medium  was  considered  to  be  isotropic. 

The  following  sections  present  a  brief  exposition  of  the  method  of 
self-similar  potentials  as  applied  to  solve  dynamic  problems  in  homogeneous, 
linearly  elastic,  anisotropic  bodies. 

2.3.1  The  Solution  of  the  Equations  of  Motion 

The  general  solution  of  the  partial  differential  equations  of  motion 
(2.15)  for  an  initial  disturbance  concentrated  at  the  origin  of  the  y>_0 
half  space  can  be  written  in  the  form 

6k  =  t-xe<-yAK(ei<_)  =  0  (k  =  1,  2)  (2.16) 

which  defines  the  characteristic  surfaces  through  the  origin  in  terms  of 
a  parameter  e  .  Here  *K(&K)  is  an  analytic  function  to  be  determined 
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later.  This  form  is  chosen  because  it  yields  a  simple  expression  for 
on  the  y=0  plane  where  the  boundary  conditions  must  be  satisfied.  On  this 
plane  e^=t/x,  i.e.  the  parameter  is  uniquely  determined  for  every  point 
and  is  the  same  for  each  solution  of  the  wave  equations.  A  more  detailed 
discussion  of  Eq.  (2.16)  is  given  in  Chapter  4  in  conjunction  with  the 
mapping  of  the  x,  y,  t  space  into  the  complex  plane. 

The  parametrization  of  the  characteristic  surfaces  by  the  quantity 
implies  that  any  analytical  functions  of  6^  must  satisfy  the  equations 
of  motion.  Let  <j>  (e  )  and  v  (e  )  be  two  functions  with  continuous  first, 

K  K  K  K 

second  and  third  derivatives.  Using  the  expressions  provided  in 
Appendix  A  and  substituting  $  =  <i>^ ( 6^ )  and  ip  =  ^(e^)  in  Eq.  (2.15),  one 
arrives  at  the  equivalent  system 

-eK[ae^  +  (c+d)  x2(0k)  -  l]  ^(e^)  - 

-xk(0k)  [(a-c)e2  +  dx^(©K)  -  1]  ^(e^)  =  0 

(2.17) 

-\(eK)  [(c+d)e2  +  bx^(eK)  -  1]  $'(eK)  + 

+eK  [de2  +  (b-c)  x2(eK)  -  1]  ^(eK)  =  0 


where  <t>'  (e  )  and  y'(e  )  denote  first  derivatives  with  respect  to  e  .  The 

fC  K  K  K  K 

condition  for  the  system  to  possess  a  non-vanishing  solution  is 


F(W>  2  bdiK  (9K>  -  (btd-UK  )X^  (6k’  + 
+  (1-ae2  )  (1-de2  )  =  0 

K  K 


(2.18) 


a  (e  )  as 

k  tc 


The  solution  of  this  biquadratic  equation  defines  the  function 
b+d  -  Le2  +  (-1)*  /Q(eJ  I*5 


a  (e  ) 

K  K 


=  + 


2bd 


(2.19) 


where 


Q(e  )  =  (b+d-Le2)2  -  4bd(l-ae2)  (1-de2)  = 


(2.20) 


=  (b-d)2  -  2M02  +  Kfoe* 


The  complete  representation  of  A  (e^)  on  a  four-sheeted  Riemann  surface 
can  be  found  in  Appendix  B.  There  exist  four  different  types  of  Riemann 
surfaces  denoted  by  I,  II,  III  and  IV  according  to  the  location  of  the 
branch  points  on  the  complex  domain.  Once  the  sign  of  the  outside  radical 
is  fixed  (see  Chapter  4),  the  two  distinct  values  of  Ak(6k)  define  by  Eq. 
(2.16)  the  variable  0^  as  the  parameter  of  the  characteristic  surfaces. 

Under  these  conditions  the  real  and  imaginary  parts  of  any  analytic  functions 
$K(eK)  and  4^(0^)  give  solutions  to  the  equations  of  motion.  If  in  a 
given  domain  the  as  function  of  x,  y  and  t  assumes  real  values,  then  any 
arbitrary  real  functions  of  6^  with  continuous  derivatives  up  to  the  third 
order  will  also  be  solutions  of  Eq.  (2.15) 

The  importance  of  the  method  of  self-similar  potentials  lies  in  the 
simplicity  of  the  determination  of  the  analytic  functions  ( 0K )  and 
4^(0^).  In  the  class  of  problems  undertaken  by  this  study,  four  of  these 
functions  are  required  to  determine  the  solution.  However,  the  homogeneous 
system  (2.17)  provides  a  value  for  the  quotient 


21 


As  a  result,  the  four  functions  <^(0^),  4^(6^),  anc*  are 

not  independent  but  for  each  value  of  the  index  k  only  one  of  the  two 
functions  4^(6^),  ^(e  )  needs  to  be  determined;  the  other  can  imnediately 
be  obtained  by  use  of  Eq.  (2.17).  The  complex  solutions  of  the 
Smirnov-Sobolev  method  must  contain  both  of  these  functions  for  each  of 


the  two  wave  fields  since  4^(0^)  is  associated  with  the  longitudinal 
component  and  ^(e^)  is  associated  with  the  transverse  component  of  the 
displacement  vector.  It  turns  out,  however,  to  be  simpler  to  pursue  the 
solution  in  terms  of  new  functions  ft  (e  )  hereafter  referred  to  as 

K  K 

"potentials."  Either  of  the  two  relations  in  Eq.  (2.17)  can  be  used 
to  define  these  potentials;  the  first  of  the  two  has  been  chosen.  Thus, 


4>'(e  ) 

_ K  K _ 

a  (e)  C(a-c)  el  +  d\zJe)  -  1] 


K  K  1C 


y'  (0  ) 

- 5 - - =  a  (e  )  (*=1,2) 

eK  [ae*  +  (c+d)  a^)  -  1]  K  K 


(2.21) 


The  ratio  of  the  functions  <t> 1  ( 0  )  and  y'(e  )  is  a  measure  of  how 

IC  K  K  K 

close  the  solution  comes  to  a  pure  longitudinal  wave  in  the  anisotropic 
media.  In  the  case  of  isotropic  materials  (Appendix  G)  the  denominator 
of  the  second  fraction  in  Eq.  (2.21)  becomes  zero  for  K=1  while  the 
denominator  of  the  first  fraction  vanishes  for  K=2.  The  longitudinal  wave 
is,  then,  described  by  4>j(e^)  and  the  transverse  wave  by  y^^) 

^(0^  =  *£(e2)  =  0. 
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2.3.2  The  Determination  of  the  Self-Similar  Potentials 

For  given  initial  or  boundary  conditions,  the  determination  of  the 
solution  of  a  particular  wave  propagation  problem  requires  the  calculation 
of  the  potentials  ft^e^).  The  methods  of  complex  analysis  are  extremely 
useful  at  this  point.  The  real  boundary  tractions,  given  as  homogeneous 
functions  of  space  and  time,  can  be  expressed  as  analytic  functions  in  the 
half-space  y>0  by  means  of  the  Schwartz  integral  formula  [9]: 


oo 


—  00 


★ 

where  T  indicates  the  real  traction  and  T  the  complex  traction  as  a 
function  of  the  complex  parameter  e  which  assumes  real  values  t/x  on  the 
boundary  y=0  plane.  The  final  result  is  given  in  terms  of  complex  functions 
where  the  real  and  imaginary  parts  correspond  respectively  to  the  real  and 


•  ,  _  ★ 
imaginary  parts  of  the  complex  boundary  conditions  given  by  T  . 

It  can  easily  be  verified  that  the  complex  valued  displacements  and 

stresses  (indicated  by  an  asterisk)  can  be  written  as  follows 

ux  =  H'f  Z1  [•2-**S(e>]  a^e)  de 

[2  Qk  m 

z  f[e2  +  x2( e)]  [ae2  +  dx2(0)  -  i]  ^(0)  de 

tc—  1  . 

0 


*/  -  ^  1 
vp  ■ 

-  2 

z 

.<=1 

(e)L  (e)  nje)  de] 

at  | 

*0  J 

*  a2  1 

-  2 

eK 

r  K 

y  at  1 

z 

.<=1 

h(e)  SK(e)  nje)  de 
0 

*  _  a2 

Vp  ■  ^7 

CM  W  II 

eK 

J  eM(c(e)  nK(e)  de] 

0  * 

(2.23) 
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where 

Lje)  =  -  [e2+A2(e)]  [d[ae2+dA2(e)  -  1]  +  c[l-dA2(e)]j 

SK(e)  =  [e2+A2(e)]  [b[ae2+dx2(e)  -  l]  -  c(c-d)e2  ]  (2.24) 

MK(e)  =  d[e2+A2(e)]  [ae2+(d-c)A2(&)  -  1] 


The  self-similarity  of  the  potentials  implies  that  the  boundary 

tractions  are  homogeneous  functions  of  the  space  and  time  variables  of 

order  n=-2.  The  generalization  of  the  method  to  obtain  solutions  for  a 

problem  where  the  given  tractions  are  homogeneous  functions  of  order  n  is 

only  a  matter  of  a  simple  n-fold  differentiation  (for  n  positive)  or 

n-fold  integration  (for  n  negative)  with  respect  to  t  (see  reference  [47], 

section  2.8).  On  the  y=0  boundary  plane  only  the  tractions  a  and  t 

y  *y 

are  applied.  The  complex  self-similar  tractions  on  the  same  plane  are 
given  by  a  double  integration  with  respect  to  t  of  the  corresponding 


This  system  of  equations  is  solved,  after  the  integral  is  removed  by  a 
differentiation  with  respect  to  e,  to  give  the  expressions  for  the 
potentials  in  terms  of  the  self-similar  tractions 


ft-|  (e) 


^(e)  = 


eM2(e)  (Ey  /p)  -  x2(a)  S2(e)  (Txy  /p) 
n(e) 

-9M-| (e)  (Ey'/p)  +  x^e)  Sl(e)  (Ty*'/p) 
n(e) 


(2.26) 


where 


n(e)  =  0  [a] (0)  S^o)  f12(0)  -  a2(@)  s2(e)  m7(o)]  = 


=  cdo  [x1 (©)-x2(0)]  [o2+x12(o)3  [0“+a2^(8)]  R(0) 


2  2, 


(2.27) 


R( 0 )  =  (K^ 02-b)  a1(0)a2(0)  +  a02-1 


(2.28) 


Here  R(0)  is  the  equivalent  of  the  Rayleigh  function  for  anisotropic  media. 
2.4  The  Anti  plane  Problem 

A  general  definition  and  the  treatment  of  a  variety  of  antiplane 
problems  can  be  found  in  reference  [23].  In  this  study,  a  special  case  of 
antiplane  problems  is  considered,  where  all  the  non-vanishing  displacement 
and  stress  components  are  parallel  to  the  normal  to  the  plane  (x,  y)  and 
are  functions  of  x,  y,  and  t  [18].  In  such  a  case,  points  lying  in  planes 
parallel  to  (x,  y)  before  deformation  cease  to  do  so  after  deformation. 

The  development  of  solutions  for  transversely  isotropic  materials,  presented 
in  the  next  chapter,  requires  the  formal  solution  to  this  kind  of  problem. 
The  equation  of  motion  is  written 

a2u  a2u  32u 

2~  +  d  2~  =  2~ 

W  dy  3t 


e 


(2.29) 


where  the  elastic  parameters  e  and  d  are  defined  from  the  stress-strain 
relations 


The  displacement  u2  is  only  a  function  of  x,  y  and  t.  It  is  not 
difficult  to  show  that  the  characteristic  surfaces  in  this  case  are 
defined  by  the  equation 

63  =  t-x03-yx3(03)  =  0  (2.31) 

where 

a3(03)  =  /i7d  (e-1  -  e| )*5  (2.32) 


The  solution  is  given  in  terms  of  a  displacement  function 
w  =  Re  [w  (03)J 

and  in  the  case  of  self-similar  tractions  is  expressed  by  the  relations 


/  3 

0w  (0)  d0  (2.33) 

0 

fyz*/p  =  -  d  w *  (0)  do 

0 

For  a  given  boundary  traction  the  third  of  Eqs.  (2.33)  determines  the 
displacement  function  w  (0)  on  the  boundary.  The  complete  solution  is 


given  by  Eq.  (2.33). 
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2. 5  Wave  Propagation  in  a  Half-Space  Resulting  from  the  Sudden  Application 
of  a  Line  Force  on  the  Surface 

The  determination  of  the  potentials  for  two-dimensional  problems  in 
anisotropic  materials  follows  the  same  sequence  of  steps  used  to  determine 
the  potentials  in  isotropic  materials  [47,  Chapter  2].  Only  a  brief 
outline  will  be  given  here.  The  magnitudes  of  the  vertical  and  horizontal 
forces  are  denoted  by  Fy  and  F^. 


2.5.1  The  Vertical  Surface  Force 

The  complex  tractions  obtained  from  the  boundary  conditions  can  be 
written  in  the  form 


*i 


r  {e)  =  -  ~ 
y  nr 


Txy*'(e)  =  0 


(2.34) 


which  determine,  by  means  of  Eq.  (2.26)  the  potentials 

Fv  [ae2  +  (d-c)  x2(e)  -  1] 

n,(e)  =  -  - K — p - 

1  l7rpC  [x^e)  -  x2(e)]  [e2+x12(e)]  R(e) 

F  [ae2  +  (d-c)  x?(e)  -  1] 

M0)  =  t-V - T - 5 - 

11Tpc  [x-j (e)  -  x2(e)]  [e2  +  x2(e)]  R(e) 


(2.35) 


The  solution  for  a  step  vertical  line  force  is  then  given  by  the  following 
expressions 
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•  Fv  ?  f  f  eVe) 


ex  (e)  [ae2+(d-c)  x?  <e)  -  1] 

^  j"  K 


wa  (e)  -  x,  (0)]  R(e) 

Q  K  OK 


de 


] 


0 

F.  2  [■  ^(l-ae2)  [(d-c)e2  +  bx^_K(e)  -  1]  _  ] 

^  [  J  [xk(0)  -  x3_|c(e)]  R(e)  de  J 


U.  =  T 


y  i,eb  K,, 


a  =  - 


F..  2  r  XK(eK)  [K-|62(l-ae2)  +  c[ae2+bx2_K(eK) 


x  i  irb 


2 

z 

c=l  L 


- 1]] 


K  C‘«<6k)  -  x3-k(Vj  R(9k> 


.  d  x  (e  )  n-bx|.K(eK)]  ] 

*  F»  2  rv^1  (,'ae?)  [bi3-K(8K>  ■ 

(c-d)2  e  2  xK(eK)  x|  (e  )  1 

+  W?  j 

*  f,  2  r  oAk-'1  <Kie<-b-c+d)  l 

vwil*;  [».(•«)- *3- A»  AH 


(2.36) 


where 


x-yx1 (e  ) 

<  K 


(2.37) 


and  x'  is  the  derivative  of  the  function  x  (e  )  with  respect  to  0 

K  K  K  K 

(see  Eq.  (4.3)). 


2.5.2  The  Horizontal  Surface  Force 

In  a  similar  manner,  the  complex  boundary  tractions  for  the  horizontal 


surface  force  are 
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Ey*'(e)  =  o 

*.  F 

T  (e)  =  -A 
xy  v  '  nr 


(2.38) 


which  determine  the  potentials 

F.  x9(e)  [b{ae2+dx|( e)  -  1}  -  c(c-d)  e2] 
Me)  =  t— 7  — - ^ - 5 — r - 

1  l7rpcd  e  [x-,  (e)  -  x2(e)J  [e2+x2(e)]  R(e) 

F  x,(e)  [b{ae2+dx2(e)  -  1}  -  c(c-d)  e2] 
M0)  =  -  7-^7-r  — - - - 5—5 - 

2  1irpcd  e  [x1(e)-x2(e)]  [e2+x2(e)]  R(e) 


(2.39) 


The  displacement  and  stress  field  for  a  horizontal  line  force  in  the  form 
or  a  step  function  in  time  applied  at  the  origin  is  given  by  the  following 


expressions 

*  Fh  2  r  Am®1  r(e)  [M“e2td4  J0)  -  1)  -  c(c-d)e2] 

ux  =  T^d  l  - — •  [>;(»)  -  >3 .pi »(.) - ^ 

.  F.  2  f  rV  (e)  [ae2+(d-c)  x2(e)  -  1]  1 

V'T^l  l  l-V(e)  -  Ve(e)j  R(e)  de  J 


.  i  2  re.  v(eK>  weK> 

°x  -  -  f,  L  K  t\A>  -  WVJ  r<®.) 


(2.40) 


*-  rA  ?  re.  v(eK>  x3-K<eK>  M  -  b-ctd)i 

°y  i-»  L  K  LMV-WVU  J 

*  fh  2  r*3-K(eK>  »C"I>2  K  [1'b>X(8X)] 

V  ■ '  i»  K.,  I  «;  iw  -  H-SV)  R<V 


m 


2.6  Reflection  of  Waves  from  a  Plane  Boundar 


When  a  disturbance  generated  by  a  source  in  the  interior  of  a  half¬ 


space  reaches  the  free  surface  the  conditions  imposed  by  the  presence  of 


the  boundary  must  be  fulfilled.  The  result  is  the  creation  of  reflected 


and  surface  waves.  Prior  to  the  arrival  of  the  incident  waves,  the  motion 


of  the  body  is  not  affected  by  the  presence  of  the  free  surface;  i.e.  the 


displacements  and  stresses  are  identical  with  the  solution  in  an  unbounded 


medium.  Upon  arrival  of  each  incident  wave  at  the  surface,  two  reflected 


waves  appear  and  their  effect  on  the  solution  is  represented  by 


corresponding  potentials  which  must  be  determined. 


Since  the  characteristic  surfaces  of  the  quasi-longitudinal  and 


quasi -transverse  incident  waves  must  pass  through  point  x=0,  y=y0.  where 
the  force  is  applied,  they  can  be  expressed  by  the  equations 


6k  5  t-x0K '  (y-y0>  =  0  (lc=1’  2) 


(2.41) 


In  the  Smirnov- Sobolev  procedure  it  is  essential  to  parametrize  the 


characteristic  surfaces  of  the  reflected  waves  in  a  way  that  each  incident 


and  its  two  reflected  waves  are  given  by  the  same  value  of  the  parameter  on 


the  free  surface  y=0.  Thus,  if  <=1  in  Eqs.  (2.41)  denotes  the  quasi¬ 


longitudinal  wave,  the  characteristic  surfaces  for  the  reflected  quasi¬ 


longitudinal  wave  are 


+  (y+yQy  *-j  (0ii )  -  0 


(2.42) 


and  the  surfaces  for  the  reflected  quasi -transverse  are 


$12  =  t-xeig  Y^2(012^  yo^(0i2)  =  ® 


(2.43) 
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Similarly,  the  reflected  quasi-longitudinal  and  quasi -transverse  waves 
produced  by  an  incident  quasi-transverse  v/ave  are  expressed  respectively 
by 

<$21  =  t-x02]  yxi  ( e2i )  +  ^0^2 ^ e21  ^  =  ^  (2.44) 

^22  =  b-x022  +  (y+y0)  ^2^22^  =  ^  (2.45) 


(2.48) 


The  double- indexed  reflected  potentials  are  related  to  the  incident 
potentials  having  the  same  first  index  while  the  second  index  indicates 
the  type  of  wave,  i.e.  <=1  refers  to  quasi-longitudinal  and  k=2  to 
quasi-transverse  waves.  For  example,  &12(0)  indicates  the  potential  of 
the  reflected  quasi-transverse  wave  caused  from  an  incident  quasi¬ 
longitudinal  wave.  The  displacements  and  stresses  for  self-similar 


potential  problems  can  be  written 


<■*  ,  J 

K-l  «-  0 


f  fK 

i  l  I  cex* 


(e)  [e2  +  x^(e)]  aK(e)  de  + 


u 

2  r  1 

e  /  cex-(e)  [e2+x?(e)]  a  .(e)  de 
j=l  J  J 


uy  ‘  '  Jt  J 

k=1  L  -'g 


!  [  1  ^2+a<(0^  [a02+d^(e)  -  1]  nK(e)  de 


■ll 


1 

[02+x2(e)]  [ae2+dx2(e)  -  1]  aKj.(e)  de  I 


*  a2  2  r  r  k 

ax/p  =  ~1  E  /  ak(0)  lk(0)  nK(0)  de  + 
3t  k=1  L 

e  . 

2  r  1 


d  r  - 

+  e  j  x j ( e)  L j ( e )  a  .(e)  de 
J=1  Jo 


*  a2  2  r  r K 

J P  =  — 2  Z  I  M0)  ^^(0)  Q,,(0)  de  + 
y  ar  <=i  L  J  K  K  K 
o 
e 

f  X .(e)  s.(e)  a  .(e)  de  1 


2 

+  E 
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*,  _  3 

Vp ' 


j,  ( / 


eMK(e)  nK(e)  de 


Z  r  KJ 

i  /  eM.(e)  n  .-(e)  de 

j=i  K  3  K3 


The  above  equations  simplify  considerably  for  points  on  the  surface,  where 


el  =  =  e12 


and 


e2  =  e21  =  e22 


(2.50) 


It  is  not  difficult  to  show  that  for  a  force  parallel  to  the  y-axis  the 
incident  potentials  are 
F. 


nf(e)  =  - 


2i^pbd  x, (e)  Cxf(e)  -  x|(e)]  [e2+A^)] 


aj(e)  =  -  v 

^  04 


1 


2i-nPbd  x2(e)  [x|(e)  -  x|( e)]  [e2+x|(e)] 


(2.51) 


and  for  a  force  parallel  to  the  x-axis 

,h  -  •  F 


n'j(e)  = 


ae2+dx2(e)  -  1 


2i7rpcd  e  [x2(e )  -  x|(e)]  (1-ae2)  [e2+x2(e)] 


nj(0)  =  ~  h 


ae2+dx2(e)  -  1 


2i7rpcd  e[x2(e)  -  x|(e)]  (1-ae2)  [e2+x2(e)] 


(2.52) 


Therefore,  the  general  solution  to  a  line  force  applied  at  x=0,  y=0  parallel 
to  the  positive  x  or  y  axis  and  varying  as  a  step  function  in  time  can  be 
determined  from  Eqs.  (2.49)  by  integrating  once  with  respect  to  time.  In 
this  way,  all  the  functions  involved  are  determined  completely. 
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2.7  Surface  Displacements  and  Stresses 

The  examination  of  the  dynamic  response  of  a  point  located  on  the 

free  surface  y=0,  due  to  a  disturbance  in  the  interior  of  a  half  space  is 

greatly  simplified  by  the  relations  in  Eq.  (2.50).  After  some  algebraic 

manipulations,  the  surface  displacements  and  stresses  for  a  step  force 

applied  at  y=yn  and  parallel  to  the  surface  can  be  written  as 
U  (H) 

u*  <h)  _  Fh  \  r  f6«  y(9>  >3-^(e)  V _ de 

x°  t,cdb  ^  l  Jo  [xic(o)-a3_k(0)]  C**(e)-x|_ic(e>3  R(e) 

umh)  .  A_  f  r  A  uyAa> _ de 

y0  mpdb  k=1  [  J  [xic(e)-A3_ic(e)]  [x^(e)-x|_ic(e)]  R(e) 


*  (H)  =  Fh 
xO  ~  iirbd 


H- 

<=ll  S' 


0.  A.  (e«)  “l"*  (9e> 

>3  R<6k> 


(2.53) 


=0 


where 

U^(e)  =  (K^-b)  [de2+bx2(e)-l]- 

-b[ae2+dx|_(c(e)-l  ]  +  ce2(K102-b-c+d) 

l/H^(e)  =  (K,e2-b-c+d)  [ae2+dx~  (e)-l]  + 
y<  I 

+c(Kie2-b)  x2_K(e)  -  c(ae2-l ) 


(2.54) 


(2.55) 


£ 

f 
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(2.60) 


The  stresses  are 


;r- 

iirdt)  K=ll  6* 


"S’  <»*> 


[x  (e  )-x~  (e  )] 


x3-,c<e*> 


[x!(e«>-xL(e.>]  R(eK> 


*(V)  *(V) 

yO  xy  u 


where 


U**v)(9e)  -  bfd  [aeR+dJ^(eK)  -  1]  +  c  [l-dV^njx 

x  [(K,e^-b)  xf(9t)  -  ao^+lj  * 

+  e2(K,e2-b-c+d)  [  cd(ae2-l )-ad  [de2+bx~  (e  ) 

K  I  K  I  <  K  J~K  K 


(2.61) 


(2.62) 


V.t*v*v<  -C<:~ 
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3.  THE  POINT  FORCE  IN  A  TRANSVERSELY 
ISOTROPIC  HALF  SPACE 


3. 1  General  Remarks 

Anisotropic  materials  possessing  an  axis  of  elastic  symmetry  at 
each  point  such  that  all  directions  perpendicular  to  that  axis  are 
equivalent  are  called  transversely  isotropic.  The  propagation  of  transient 
disturbances  in  media  of  this  type  is  a  subject  of  both  mathematical  and 
geophysical  interest.  Highly  compressed  sedimentary  material  where  the 
vertical  compression  induces  significantly  different  properties  in  that 
direction  can  be  modeled  as  transversely  isotropic.  Of  particular 
importance  in  seismology  and  earthquake  engineering  are  composite  media 
having  different  isotropic  layers.  When  the  thickness  of  individual 
layers  is  small  compared  with  wave-length,  many  of  the  characteristics 
of  elastic  wave  propagation  are  similar  to  those  of  transversely  isotropic 
media.  Backus  [4]  has  given  conditions  under  which  the  homogeneous, 
transversely  isotropic  medium  can  indeed  be  the  long-wave  equivalent  of  a 
layered  isotropic  medium  made  of  just  two  homogeneous,  isotropic  materials. 

The  high  degree  of  symmetry  exhibited  by  transversely  isotropic 
materials  reduces  the  number  of  independent  elastic  constants  required 
to  define  its  behavior  to  five  [21, §105].  The  same  number  of  constants  is 
also  required  by  hexagonal  crystals  which  have  an  axis  of  elastic  symmetry 
of  order  six.  Hence,  "transverse  isotropy"  and  "hexagonal  aeolotropy" 
are  used  interchangeably. 

Very  few  solutions  for  any  but  trivial  problems  are  available  in  the 
literature  even  for  isotropic  materials  in  a  three-dimensional  space.  This 
study  utilitizes  the  method  of  rotational  superposition  of  plane  and  anti- 
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plane  problems  to  obtain  solutions  for  the  problem  of  a  suddenly  applied 
point  force  on  the  surface  of  a  transversely  isotropic  half-space.  The 
axis  of  elastic  symmetry  is  taken  perpendicular  to  the  free  surface.  The 
case  of  the  subsurface  load,  although  not  treated  here,  can  be  obtained 
by  utilizing  the  same  technique  on  the  two-dimensional  solutions  for  the 
point  force  in  the  interior  of  the  half-space  given  at  the  end  of 
Chapter  2. 

3.2  Equations  of  Motion 

Consider  a  cylindrical  coordinate  system  (r,  u>,  y)  where  the  y-axis 
is  parallel  to  the  axis  of  elastic  symmetry.  The  stress-strain  relations 
for  a  transversely  isotropic  material  can  be  written  as 


°r  =  C11  er  +  c13  £o>  +  C1 2  ey 

°a)  =  c13  er  +  C11  Gu)  +  C1 2  Ey 

ay  =  c12  (er+eu>}  +  c22  ey 

Vy  =  c66  Yry 

C66  Yyto 


The  equations  of  equilibrium  for  the  medium  are 


3o„  ,  3t  3t  o-o 

-X  +  1  — Ea  +  -JOL  +  JLJS.  +  z  =0 
3r  r  3a)  3y  r  r 


(3.1) 


m 
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In  addition  to  the  elastic  parameters  defined  by  Eq.  (2.7)  a  further 


constant  "e"  is  often  introduced, 
cll‘c13 

e  ‘  2p 


(3.5) 


Solutions  will  be  given  in  terms  of  the  five  elastic  parameters  a,  b,  c, 
d  and  e  which  must  satisfy  the  conditions 

a>0,  b>0,  d>0,  e>0,  b(a-e)  -  (c-d)2>0  (3.6) 


These  are  equivalent  to  the  requirement  that  the  strain  energy  be 
positive. 


3.3  The  Method  of  Rotational  Superposition 

The  application  of  the  method  of  rotational  superposition  to  solve 
certain  three-dimensional  wave  propagation  problems  in  isotropic  materials 
can  be  found  in  references  [12,  18,  47].  Some  repetition  of  the  basic 
principles  cannot  be  avoided  in  applying  the  method  to  anisotropic  materials. 

Consider  the  half  space  of  Fig.  1,  where  (X,  Y)  defines  the  plane 
of  deformation  for  the  plane  strain  and  antiplane  problems,  and  (x,  y)  is 
the  reference  plane  for  the  coordinate  angle  w.  If  U^,  Uy  denote  the 
displacements  of  the  plane  problem  and  the  displacement  of  the  anti  plane 
problem,  their  contribution  to  the  three  dimensional  displacements  at 
point  P  (r,oa,y)  can  be  expressed  by  the  relations 

Ur  =  Uy  [rCOS  (u>q-u>),  Y,  t]  cos  (u>q-u>)  - 
-  Uz  [rcos  (ujq-u>),  Y,  t]  sin  (u>q-u>) 

Ujj  =  [rcos  U0-«),  Y,  t]  sin  (u>q-o>)  + 
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+  Uz  [rcos  (oig-o)),  Y,  t]  cos  (aig-w)  (3.7) 

uy  “  UY  [rcos  (a)0-u),  Y,  t] 


The  superposition  of  the  effects  of  all  plane  and  antiplane  problems  at 
different  angular  orientations  Wg  multiplied  by  weighting  functions  ^(“g) 

and  f2(o)g)  which  are  going  to  be  defined  later,  results  in  a  three- 
dimensional  wave  field  for  Isotropic  materials  [18].  The  same  principle 
applied  in  the  present  case  gives  the  expressions 


u 


r 


u 

0) 


[rcos  ( cjg— (i) ) ,  Y,  t]  cos  (ti)g— (i) )  f | 

[rcos  (oig-o)),  Y,  t]  sin  (a>g-a>)  f2 

[rcos  (ojg-oj),  Y,  t]  sin  (tug-w)  f-j 

[rcos  ((Dg-oi),  Y,  t]  cos  (wg-w)  f2 

[rcos  (oig-u),  Y,  t]  fj  (ug)  dug 


(“g)  du)0  ' 

(<i>g)  dujg 

Ug)  d(Dg  + 

(Oig)  dWg 


(3.8) 


Straightforward  application  of  these  displacements  to  the  equations  of 
motion  (3.4)  shows  that  the  latter  are  satisfied  if 


32U, 


'11 


+  c 


3X 

(C12+<W 


66 


2  2  2 
d\  a  UY  3  UY 

+  ^c12+c66^  ~  =  P  ~ 


3Y‘ 


3X3Y 


3t 


3X3Y 


32U, 


+  C 


66 


3X 


32Uy  32Uy 

+  C22  TX  =  0  TX 


3Y 


3t 


(3.9a) 


and 

cn~ci3  *2uz  _  &2uz 

2  3X2  66  3Y2  "  at2 


(3.9b) 


1 
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These  are  the  equations  of  motion  for  the  plane  strain  and  antiplane 
problems.  Consequently,  the  rotational  superposition  of  plane  solutions 
results  in  a  three-dimensional  field  which  satisfies  the  equations  of 
motion  and  thus,  by  uniqueness,  is  the  solution.  It  is  interesting 
to  see  what  happens  to  stresses  under  this  operation.  By  use  of  Eqs.  (3.1), 
(3.3)  and  (3.8)  it  is  not  difficult  to  arrive  at  the  expressions 


,1T. 


ft  2  2  UWX 

or  —  I  |I“^ii  cos  (ojq-w)  +  c^ ^  sin  (w^-w)l  a y 

0 

3UV 


9U) 

3X 


ouy  \ 

c12  TT  |  fl  ^“0^  dw0  ' 

(cn‘ci ?)  f”  3U7 

- j— —  J  -gjr  sin  2(wq-o))  f2 (^q )  dcJo 


r,  2  2  BUX 

=  J  ([cll  Sin  (“0"o,)  +  C13  C0S  (w0_a,)]  1T  + 


3Uy 

+  ci2  ir!  W  daio + 

7T 

Cn-C,,  r  3U, 

J  13)  Z  „  • . 


+  - 2 -  !  lx”  S3n  ^“0^  du0 

0 

rv  3UX  3Uy 

=  J  *-C12  TT+  c22  TH  M'V  dw0 

0 

-j  cos  (uq-u)  ^(wg)  d 


f  1 

r  3  u  w  sUw 

'  J  c 66 

1 _ p.  +  _ L 

¥  3X 

‘ 

311 

J  C66 

IT  si"  <■ 

0 

-  fc 

'3UX  3Uy 

"  J  C66 

IT  +  lx" 

'2V  O' 


f  3UZ 

+  JC66  IT  cos(b,0_0,)  f2^0)  du0 


(3.10) 


7T 

f  cn-ci3 

Tra)  =  l  -  ~  sin  2(a,0-a))  f^g)  du>g  + 

IT 

f  ^^7 

+  J  ~  2  ~  IT  cos  2(V“}  f2(“0)  S 

0 

By  virtue  of  the  relations  between  stresses,  strains  and  displace¬ 
ments  in  the  two-dimensional  case,  these  equations  can  be  rewritten 

f 71  •  2 

=  Ey  f-|  (wq)  doTg  -  ( Cl  1  -  C-]  3  )  j  Ey  sin  (u)g-Clj)  fl(ojg)  d(l)g  - 

(,1T  0 

TyZ  sin  2(wg-oi)  f 2 ^ ^g 5  d“g 

Y  *  f* 

°u>  =  Jh  f1^0^  da30  "  ^Cll"C13^  I  h  COs2(uq-u)  f-|(wg)  dajg  + 

+  J  Tyz  sin  2((Ug-aj)  f 2 ( ^g )  da)g 
0 

IT 

Oy  -  JE y  f-j  (<»)g)  daig  (3. 

0  7T 

Try  =  /  TXY  fl((i)0}  dlli0  " 

V 

"  YZ  s1'n(Vw)  ^2^U(P  d(i)0 

0  7T 

V  =  TXY  Sin(a)g-O))  ^((Og)  d0)g  + 

'  -TT 

+  TyZ  cos(ojg-a))  f2(“g)  d“g 

c°-c  r" 

Tru>  =  ~ 2  13  J  EX  S1'n  2(V°j)  fl(a)0}  d“0  + 

_  7T  0 

J  TyZ  2 (aig-u)  ^2^0^  duig 


(3.11) 


Here,  Ey,  Ey,  Tyy,  Ty^,  Tyz  denote  the  stresses  and  Ey  the  strain  in  the 

X-direction  of  the  two-dimensional  problems.  It  can  be  seen  that  ay, 

t  and  x  are  formed  by  a  rotational  superoosition  of  the  two-dimensional 
ry  yu 
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stresses,  while  the  remaining  stresses  involve  additional  terms.  This  is 
a  very  important  observation  because  the  three  quantities  o  ,  x  ,  x 

J  I  J  J  W 

must  be  specified  on  the  boundary  to  define  the  three-dimensional  problem. 
Thus,  by  inversion  of  the  three  integral  equations  involved,  the  two- 
dimensional  quantities  can  be  determined. 

The  rotational  superposition  of  plane  strain  and  anti  plane  problems 
in  anisotropic  media  must  be  differentiated  from  the  same  method  as 
applied  in  the  case  of  isotropy.  There,  the  medium  for  each  contributing 
plane  or  anti  plane  solution  is  the  same  and  the  superposition  of  the  wave 
fields  is  a  natural  consequence.  On  the  contrary,  the  solutions  that  have 
been  superimposed  in  this  study  correspond  to  two  quite  different  materials. 
It  is  a  remarkable  fact  that  the  rotational  superposition  of  the  plane 
displacement  fields  leads  to  expressions  for  the  stresses  oy,  xry,  and  xy^ 
which  can  be  solved  easily  to  define  the  two-dimensional  problems. 

3.4  Axi synwetri c  Problems 

A  three-dimensional  axisymmetric  problem  can  be  solved  by  a  rotational 
superposition  of  the  solution  of  a  unique  plane  strain  problem  [12,18]. 

The  determination  of  the  plane  problem  is  simplified  by  observing  that 
the  two-dimensional  field  must  be  symmetric.  Here,  the  weighting  functions 
are  taken  as 


f-|  (w0)  =  1 


f2(u>0)  =  0 


(3.12) 
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After  introduction  of  the  new  variable  W  =  u  -  the  displacements  and 
stresses  are  rewritten  as 


7T 

ur  =^UX  [rcosW,  Y,  t]  cosW  dW 


u  =  0 

O) 


(3.13) 


71 

=  fuy  [rcosW,  Y,  t]  dW 
tv 

=  f (sx  [rcosW,  Y,  t]  - 
*0 

-  (c-j ^ -c^ 3 )  Ex  [rcosW,  Y,  t]  sin2w)  dW 

t r 

=  JiEx  trcosW»  Y»  *3  ‘  (cn-c13)  Ex  [rcosW,  Y,  t]cos2W ) 

7T 

=  jZy  [rcosW,  Y,t]  dW 

*4 

=  /  TXY  [rcosW,  Y,  t]  cosW  dW 


T  =  T  =0 

ru  yw 

An  especially  convenient  form  of  the  boundary  tractions  and  xry  on  the 
y=0  plane  can  be  obtained  by  using  the  transformation  u=rcosW: 


<jy(r,  0,  t)  -  2  J  £y(u,  0, 


/r2-u2 


,(r,  0,  t)  =  \  J  Txy(u,  0,  t) 


(3.14) 
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These  are  integral  equations  of  the  Abel  type  [14].  Their  solutions  are 

udo 

cyo,  u,  x;  - 

0 


yx,  o •  v  -  «yt».  o,  t) 


JF7 
2 


(3.15) 


TXV(X’  <>•*>=  IT  4  J\ ‘>7 


do 


#7 


which  give  the  boundary  tractions  of  the  plane  strain  problem  in  terms 
of  those  of  the  original  three-dimensional  problem.  The  methods  of 
Chapter  2  are  now  applicable  to  solve  the  two-dimensional  problem.  The 
final  three-dimensional  case  is  found  by  integrating  this  solution  as  in 
Eq.  (3.13). 


3.4.1  The  Vertical  Point  Force 

Consider  a  transversely  isotropic  half-space  subjected  to  a  normal 
stress  aQ  applied  over  a  circular  region  expanding  in  tine.  The  speed  of 
expansion  is  a  and  the  center  coincides  with  the  origin  of  the  coordinate 
system.  The  boundary  conditions  can  be  written 


CTy  =  -a0 


=  0 


T  =  T  =  0 
yr  yw 


for 

for 

for 


r  at 
r  >  at 

0  <  r  < 


)  0  <_u  <  2ir 

y=0 


(3.16) 


and  correspond  to  a  loading  with  a  total  vertical  force  V  =  ira2t2OQ  = 
2  2 

Fyt  where  Fy  =  ira  Oq.  By  means  of  Eqs.  (c 
defined  by  the  conditions  on  the  y=0  plane 


2  2 

Fyt  where  Fy  =  ira  <Jq.  By  means  of  Eqs.  (3.15)  the  plane  strain  problem  is 
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The  complex  displacements  of  the  two-dimensional  problem  of  a  step  force 
of  magnitude  Fy  can  be  written  from  Eqs.  (2.23)  by  differentiating  twice 
with  respect  to  time  and  multiplying  by  a  factor  of  1/2.  Thus, 


A.  ?  K+(d-o  >Ue*>  -  o 

2„Zp  V  -x3^e«u  R(9,) 

Fv  2  e,(aeZ-1)  [(d-c)i>2+l»f.|C(i>|C)  -  1] 
2i.20b  K>1  a;  Cx|c( e^)  -  a3.k(8,.)]  R(eK) 


(3.21) 


In  the  three  dimensional  problems  the  expressions  a^  and  a/, 
equivalent  to  the  two-dimensional  6  and  6*  are  written 

K  K 


a  =  t  -  rcosW  e  -  yx  (e  )  =  0 

KT  K  K  K 


(3.22) 


a'  =  -  rcosW  -  yx* (e  ) 


<  1C 


(3.23) 


The  displacement  field  can  be  obtained  from  Eqs.  (3.21)  and  (3.13).  A 
change  of  the  variable  of  integration  from  W  to  results  in  the 
following  expressions 


*  2  f  j  [ae^(J-c)  >|-c(aK)  -  0 

'  2.2p  .-I  J  V8*’ 

Ck 


de 


* 

u  =  0 

0) 


Fv  ?  r  [(d-c)  e<2+biL'9,)-i] 


^  2ir2pb  k 


/ 

C< 


Gp(0p)  E\,<0  -  W8*”  r<6k>  d9e 


where 


G  (e  )  =  re  sinW  =  [r2e2  -  [t-yx  (e  J]2]*5 


(3.24) 


K  K 


K  K 


(3.25) 


and  C  are  the  contours  of  integration  on  the  complex  plane.  The  stresses 


can  be  determined  in  a  similar  manner  but  will  not  be  given  here. 
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3.5  Ce rta in  Non- axi symmetric  Problems 


The  solutions  of  certain  dynamic  problems  which  do  not  exhibit  axial 
symmetry  can  be  formed  by  a  rotational  superposition  of  solutions  for 
both  plane  and  antiplane  problems.  The  plane  strain  problem  is  assumed 
to  be  antisymmetric  about  the  y  axis,  while  the  antiplane  problem  is 
symmetric.  Use  of  the  weighting  functions  f^  (ug)  =  coswg  and  f2(u)g)  = 
sinug,  and  the  change  of  variable  W  =  oj-oig  in  Eqs.  (3.8)  and  (3.11)  leads  to 


ur  /  cosoi  =  j cos2WdW  -  J  U^  sin2WdW 


u  /  si 

GJ 


r 

mu  =  -  | 


Ux  sin2WdW  +  /  cos2WdW 


r 

uy  /  cosu  =  J  Uy  cosWdW 

Y  f *  2 

or  /  cosu  =  J  ExcosWdW  +  (cn-c13)  i  E^sin  WcosWdW  - 


Txz  sin2WcosWdW 


o  /  cosco  =  /  zY  cosWdW  -  (c, ,-c,,)  /  EY  cosJWdW  + 


11  '"\3I  J  CX 
0 


+  J  Txz  sin  WcosWdW 


Oy  /  cos 


“  =  J 


Zy  cosWdW 


0  -Tt  —  7T 

Tyr  /  cosoi  =  J  TyX  cos2WdW  -  J  TyZ  sin2WdW 

0  ft  *0  - 

t  /  sinoi  =  -  j  Tyx  sin2WdW  +  I  TyZ  cos2WdW 

*4  _  J0 

tru  /  sinoi  =  -  (c^-c^)  /  Ex  sin2WcosWdW  + 

.  Jo 

+  j  ^XZ  cos^(cos2W-sin2W)  dW 


(3.26) 
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The  three  integral  equations  which  express  the  tractions  oy,  iyr  and  tyw 
can  be  transformed  by  using  the  new  variable  u=r  cosW  and  the  condition 

TXY  : 


TyZ  =  T  to  the  following  expressions  on  the  y=0  plane 
r 

\  I  £y(u,  0,  t) 


ay(r, 0,t) 


COSw 


0 


udu 

JJ~ 2 

/r  -u 


—  -  2  /  T(u.O.t) 

A 


du 


COSO) 


sim> 


JF? 


(3.27) 


The  relation  between  the  stresses  T^y  and  Tyz  which  produces  the  second 

of  Eqs.  (3.27)  holds  true  for  the  problems  of  this  study  (see  next 

section).  The  boundary  tractions  of  the  plane  and  antiplane  problems 

are  given  from  the  solution  of  Eq.  (3.27)  as 

X 


7TxJ 


T(X,0,t)  =  7 


rX  Tyr^u,0,t^ 


1  _a_  f  ~yu) 
"  "  TT  a  X  J 
0 


udu 

cos“ 

Tyu)(o*°>t)  udu 


(3.28) 


sirid) 


3.5.1  The  Tangential  Point  Force 

The  same  procedure  is  used  as  in  section  3.4.1.  The  horizontal 
Fh 

traction  tq  =  — ^  acts  over  an  expanded  circular  area  r  _<  at  resulting 

tra^  2 

in  a  total  horizontal  force  Fht  .  The  boundary  conditions  can  be  written 


on  the  y=0  plane  as 
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II 

o 

for 

0  1 

r  <  »' 

t 

> 

i 

yr  '  "T0COSa3  1 

i 

► 

for 

r  < 

at 

\  0  <_u)  <  2n 

ryw  =  T0sinw  J 

y=0 

yr  =  Tya,  =  0 

for 

r  > 

at 

(3.29) 


The  self-similar  tractions  of  the  plane  strain  problem  are  given  from 
Eqs.  (3.28).  Following  the  same  steps  as  in  section  3.4.1  the  solution  for 
a  step  horizontal  force  of  magnitude  can  be  written 


*  2  e  a  (e  )  [e2+A2(e  )]  a  (e  ) 

II  _  C  „  K  <'  L  K  K  <  <  K 

ux  -  -  _  z  - p 

K=1  < 


*  _  I  \  CeWteJ]  CaeWteJ  -  1] 

E  r 


UY  =  2 


where 


^(e)  - - 2 

IT 


K  =  1 

Fh6  A2(6)  S2(0) 


n(e) 


F.e  a, (e)  s,(e) 
fl2<e)  ■  -T  --Tef- 


(3.30) 


(3.31) 


It  is  not  difficult  to  show  that  the  solution  of  the  antiplane  problem  is 
given  by 

Uz  =  -  — p - j - - -  (3.32) 

2v  pd  A3(e3) 

where 

a3  =  t-e3rcosW  -  yA3(e3)  =  0 


(3.33) 


and 
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=  -  rcosW  -  yA3(e3) 


(3.34) 


Here,  *3(63)  is  given  by  Eq.  (2.32)  while 


A'(e3)  =  -  /W 


0- 


<«"  -  ^ 


(3.35) 


Direct  substitution  of  Eqs.  (3.30)  and  (3.32)  to  Eqs.  (3.26)  with  a 
change  of  the  variable  of  integration  from  W  to  yields 

2 

£ 


ur*/cos“  ■  7TT2 


2»pd  r 


2  f  AK(8K}  A3-k(0K}  [t_yAK(0K);I' 

4  4  K  W  i»A>  - 


x  [t(ae2+dA~  (0  )-l }  -  c(c-d)  e2]  de  - 

K  K  K  K 


■l 


63(63) 

e3  a3^03^  d  3 


u  /sinw  = 


w  '  2ir2pdr2 


2  r  yyyAj  w  , 

c=l4  e<  t\(eK)-A3.<(e<)j  R(bk) 


(3.36) 


x  [b{ae2+dA?  (e  )  -  1}  -  c(c-d)  e2]  de  + 

K  «3"K  K  K  K 
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[t-y  A3(e3 


C3  03  ^3^03^  A3 


>]2  1 
^rde3 


*,  Fh  ?  f  8<  J3 -A>  Caef-Hd-c)  A2(eK)  -  1] 

u  /coso)  -  — o  2  /  ?  m  rr  / „  '\  ■  n  57 «“i  x 

y  2ir  pr  <=1JC> 


'c<  W  ^teK)-A3.Kl0K)j  R(eK) 


x  Ct-yA_(e_)]  de 


where 


g3(63)  =  e3rsinW  =  [r2e2  -[t-yx 3(e3)]2]ls 


(3.37) 


The  contours  of  integration  ,  C2  and  C3  are  traced  on  the  complex  plane 
in  a  manner  to  be  described  in  Chapter  5. 


4.  THE  REPRESENTATION  OF  THE  WAVE  FIELDS  ONTO  THE 
COMPLEX  PLANES.  GEOMETRY  OF  THE  WAVE  FRONTS 

4.1  General  Remarks 

The  numerical  computation  of  the  dynamic  response  of  an  anisotropic 
medium  caused  by  an  applied  step  force  requires  the  evaluation  of  the 
real  part  of  the  complex  displacements  and  stresses  given  in  the  preceding 
two  chapters.  In  the  process,  the  need  arises  for  an  accurate  represen¬ 
tation  of  the  physical  space  on  the  complex  domain  composed  of  a  four- 
sheeted  Riemann  surface  (Appendix  B).  The  present  chapter  attempts  to 
clarify  the  regions  on  this  surface  which  represent  the  quasi-longitudinal 
and  quasi-transverse  wave  fields.  The  discussion  concerns  the  solutions 
of  the  two-dimensional  wave  propagation  problems  which  form  the  basis 
for  developing  the  transversely  isotropic  solutions. 

The  correspondence  between  points  of  the  (x,  y)  space  and  the 
complex  domain  is  given  by  Eqs.  (2.16)  which  parametrize  the  characteristic 
surfaces  of  the  solutions  to  the  wave  equations  by  means  of  the  complex 
variable  o^.  The  major  difference  from  the  isotropic  case  is  caused  by 
the  behavior  of  the  multi-valued  function  *K(0K)  which,  in  turn,  depends 
upon  certain  inequality  among  the  elastic  parameters.  In  has  been  found 
convenient  to  separate  orthotropic  materials  into  a  number  of 
different  cases,  which  greatly  facilitates  the  total  exposition. 

No  real  conceptual  difficulties  are  introduced  by  anisotropy  in 
wave  propagation  problems.  The  system  of  partial  differential  equations 
that  describe  the  motion  has  two  solutions  corresponding  to  two  types 
of  propagating  waves.  The  wave  associated  with  the  larger  velocity  is 
"nearly"  longitudinal  while  the  wave  which  travels  with  the  smaller  velocity 
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Instead,  the  curves  will  be  traced  in  terms  of  the  parameter  e  . 

Let  be  the  angle  that  the  normal  vector  n^  to  the  tangent  plane 
at  point  P  np)  makes  with  the  positive  n-axis  (Fig.  2),  and 
the  angle  that  the  ray  OP  makes  with  the  same  axis.  The  following 
relations  can  be  easily  verified 


6 


tan<V  *  TliT 

K  K 

(4.5) 

tanp  =  -x 1 (e  ) 

K  K  K 

(4.6) 

The  half  plane  n  >_  0  will  be  represented  on  the  upper  half  of  the  complex 
plane  Ime^  >_  0.  It  is  important  to  note  that  this  choice  in  the 
mapping  determines  a  corresponding  choice  of  the  sign  of  the  outside 
radical  of  the  function  A^e^).  This  is  equivalent  to  selecting  the 
two  sheets  of  the  four-sheeted  Riemann  surface  on  which  the  functions 
giving  the  solution  will  be  defined.  In  this  study,  sheets  S-j  and  S2 
are  chosen  to  represent  the  domain  of  the  complex  solutions.  The 
index  k  will  assume  the  values  1  or  2  according  to  which  sheet  of  the 
Riemann  surface  contains  the  argument  of  the  functions  under  consideration. 

Let  S-|  indicate  initially  that  complex  domain  of  the  quasi - 
longitudinal  wave  which  possesses  a  branch  cut  between  points  +  1  /»^a 
on  the  real  axis  of  9-|  (Fig.  3).  Every  point  along  the  sides  of  the 
branch  cut  is  the  image  of  all  points  lying  on  a  half  plane  tangent  to 
the  wave  front.  Thus,  the  functions  that  express  the  complex  solutions 
remain  constant  on  these  half  planes  which  are  chosen  in  such  a  way  that 
they  do  not  cross  each  other  in  each  of  the  n  >  0  or  n  <  0  half  planes. 
Consequently,  the  same  expressions  that  give  the  solution  for  every 
point  inside  the  wave  front  can  be  used  to  give  the  solution  for  points 
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outside  the  wave  front  provided  that  the  continuation  is  taken  along  the 
previously  mentioned  half  tangent  planes  [41,  §  53].  The  mapping  of  the 
space  outside  the  quasi-longitudinal  wave  front  on  the  branch  cut 
(-  1/  /a,  l//a)  is  given  in  Fig.  3.  No  point  located  in  the  region 
|  n  |  >  Jb  has  a  representation  on  the  complex  plane.  Similar 
remarks  can  be  made  for  the  quasi-transverse  wave  field  although  the 
presence  of  cusps  needs  some  extended  analysis.  Any  further  discussion 
will  concentrate  on  the  c  >  fl,  n  >_0  quarter  plane  since  the  indicated 
symmetry  makes  the  remaining  quarter  planes  equivalent. 

The  parametric  variation  of  the  angles  ^(e^)  and  0^(0^)  can  be 
examined  by  means  of  Eqs.  (4.5)  and  (4.6)  and  the  derivatives  of  the 
tangents  with  respect  to  e^;  i.e. 

(e) 


(tana  )'  =  - « - 

K  2bdA3(e  ) 


(4.7) 


K  K 


where 


and 


D  (e  ) 

(tanB  )  =  p  o  3 

K  4b Va  3(e  ) 

<  K 

(b  -  d)2  -  Me  2 
4(9)=  b  +  d+(-l)K  — - — 

K  K  /qTFT 

D  (e  )  =  <j>  (e  )♦  (e  )  +  2bde  a2  (e  )^'(e  )  = 

K  K  K  K  K  K  K  K  K  1C  K 

-  2bdt*K(6.>te^<e.>3\2  -  °2*2<eK> 


8bdc  N,e 

K  (e  )  *  (-DK  — =r 


(°K  > 

Consider  also  the  following  expressions 


(4.8) 

(4.9) 


(4.10) 

(4.11) 
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<j>'(e  )  =  - 

K  K 


(-1)K  Bbdc2^^ 


/q3  (eK) 

D1  (0  )  =  (-1)K  48b2d2c2N,0  \*(Q  ) 

K  K  1  K  K  '  K  7 


2,„  v  ^b"d)2  "  K1K28k 


4Vj 


(4.12) 


(4.13) 


An  examination  of  the  behavior  of  and  ^  in  the  interval  0  <_  ©1  ±  1  //a~ 
leads  to  the  conclusion  that  both  angles  vary  continuously  and  are 
monotonical ly  increasing  in  the  interval  (0,  7t/2)  [30 ].  Therefore,  the 
shape  of  the  quasi-longitudinal  wave  front  is  always  a  strictly  convex, 
closed  curve. 

A  complete  analysis  of  the  variation  of  the  angles  a2  and  e2 
associated  with  the  quasi -transverse  wave  is  given  in  Appendix  C.  Table  1 
presents  a  summary  of  this  investigation.  It  can  be  seen  that  the 
orthotropic  materials  under  consideration  can  be  separated  into  five 
distinct  cases,  indicated  by  the  letters  A,  B,  C,  D  and  E,  according  to 
the  existence  or  nonexistence  of  cuspidal  edges  and  the  orientation  of 
the  cusps  relative  to  the  coordinate  axes  [Table  1 ].  The  inequalities 
involving  the  elastic  parameters  associated  with  each  material  case  are 
not  always  compatible  with  the  conditions  which  accompany  each  of  the 
four  types  of  complex  planes  I,  II,  III  or  IV  (see  Appendix  B).  An 
examination  of  the  range  of  validity  of  each  inequality  results  in  the 
following  table. 


i  Uiwrithi 
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— 

Material  Case 

Compatible  Complex  Plane 

A 

I,  II,  III,  IV 

B 

I 

C 

II 

0 

II,  III 

E 

II,  IV 

In  this  chapter  and  the  next,  a  given  material  will  be  classified  by  a 
letter  indicating  the  material  case  followed  by  a  Roman  numeral 
indicating  the  type  of  the  complex  plane.  For  example,  material  A-IV 
has  elastic  parameters  that  satisfy  the  conditions  associated  with 
material  case  A  and  type  of  complex  plane  IV,  i.e.  N1  <  0,  N2  >  0,  >  0, 

N4  >  0,  K2  >  0,  M  >  0. 

4.3  The  Representation  of  the  Wave  Fields 

It  has  been  mentioned  that  the  functions  determining  the  complex 
solutions  for  the  two  types  of  propagating  waves  will  be  defined  on  the 
sheets  S-|  and  S2  of  the  Riemann  surface.  In  anisotropic  materials,  each 
sheet  does  not,  in  general,  exclusively  represent  the  field  of  a  single 
type  of  wave.  A  careful  examination  is  needed  to  define  the  areas  on  each 
sheet  that  correspond  to  the  wave  field  of  the  quasi-longi tudinal  or  the 
quasi-transverse  type  of  disturbance.  Before  proceeding  any  further  on 
this  subject,  it  is  useful  to  consider  briefly  Eq.  (4.1)  written  in  an 
equivalent  form  where  the  radicals  of  function  X^  (6^.)  have  been 


eliminated,  i.e.. 
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[bdC4  +  U2n2  +  adn4]  04  -  2?(2bd?2  +  Ln2)  e3  - 

K  tc 

-  [(a  +  d)n^  +  (b  +  d)c2n^  -  6bd£2  -  Ln^]  0^  + 

+  2?[(b  +  d)n^  -  2bd]e^  +  [n4  -  (b  +  d)n^  +  bd]  =  0  (4.14) 

This  equation,  in  general,  possesses  four  roots  in  for  every 
point  (S.n).  When  the  point  is  located  in  the  interior  of  both  wave 
fronts,  the  roots  form  two  pairs  of  complex  conjugate  numbers.  The 
value  of  for  each  pair  which  represents  the  solution  must  be  located 
in  the  appropriate  section  of  the  complex  plane  in  agreement  with  the 
mapping  convention.  The  existence  of  real  roots  for  is 
equivalent  to  saying  that  real  tangents  can  be  drawn  from  the  given 
point  to  the  wave  front.  In  accordance  with  the  method  of  the  preceding 
section  for  continuing  the  solution  along  half  tangent  planes,  the  real 
roots  must  be  chosen  appropriately.  Any  point  located  outside  both  wave 
fronts  has  four  real  roots  for  while  any  point  between  the  two  wave 
fronts  has  two  real  roots  for  e^.  Consider  now  points  located  inside  the 
triangularly  shaped  area  of  a  cusped  quasi-transverse  wave  front 
(Fig.  4).  It  can  be  seen  that  four  real  tangents  can  be  drawn  to  the 
wave  front  indicating  the  existence  of  four  real  roots  of  Eq.  (4.14). 

The  solution  is  represented  now  by  the  two  real  roots  which  are  associated 
with  the  half  tangent  planes  in  the  chosen  direction.  One  interpretation 
of  this  fact  is  that  there  are  no  roots  associated  with  the  quasi¬ 
longitudinal  wave.  In  other  words,  the  inside  of  the  cusped  region  is  a 
"hole"  in  the  quasi-longitudinal  wave. 

The  correpsondence  between  points  of  the  wave  fields  and  points 
on  sheets  and  of  the  Riemann  surface  is  expressed  by  Eq.  (4.1). 
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Let 


9  =  £  +  1? 

K  K  K 


a  (e  )  =  E  +  iZ 

K  K  K  K 


(4.15) 


then,  the  following  expressions  can  be  obtained 

Z 


e  Z  -  c  E 

K  K  K  K 


c  z  -  e  e 

K  K  K  K 


(4.16) 


These  will  be  used  in  the  next  sections  to  discuss  each  type  of  complex 
Riemann  surface  separately. 


4.3.1  Type  1  Riemann  Surface 

A  typical  configuration  of  the  complex  domain  of  Type  I  is  given 
by  two  sheets  of  the  Riemann  surface  shown  in  Fig.  5.  The  branch  cut 
(-1/t/a,  l//a)  of  the  S-j  sheet  represents  the  quasi-longitudinal  wave  front 
and  the  branch  cut  (-1//?,  l//d)  represents  the  quasi -transverse  wave  front. 
Sheets  S.j  and  S2  are  interconnected  along  branch  cuts  parallel  to  the 
imaginary  -  axis  from  points  +  eQ1  and  +  0O2  to  +  «  (Appendix  B). 

Let  e  assume  real  values  not  located  on  the  branch  cuts.  Then, 

K 

from  Eqs.  (4.15)  and  (4.16)  it  can  be  seen  that 


i.e.  as  0k  varies  from  l//a  (or  l/v/d)  to  infinity,  the  corresponding  points 
of  the  wave  field  lie  on  the  5-axis  from  the  intersection  with  the  wave 
front  5  =  /a  (or  5  =  »^)  to  the  origin.  For  purely  imaginary  values  of 
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e  (i.e.  0  =  ic  )  the  function  x  (e  )  can  be  written 


where 


'b+d+U2  +  (-1)K  v^TTFT  I  1/2 


x  (ic  )  = 

K  K 


Q(icK)  =  (b+d+U*)2  -  4bd(l+a;2)  (1+d?2  )  = 


(4.18) 


=  (b-d)2  +  2M;2  + 


(4.19) 


The  polynomial  Q  (i£K)  does  not  have  any  real  roots  in  and  as  a  result 
remains  always  real  and  positive  in  the  range  of  c  .  It  is  obvious  from 


Eq.  (4.19)  that 


Q(Hk)  <  (b+d+U2)2 


(4.20) 


Then,  from  Eq.  (4.18)  it  can  be  seen  that  XK(UK)  is  a  real  quantity. 
Equations  (4.16)  now  can  be  written 


5  =  0 


(4.21) 


It  is  not  difficult  to  show  that  X^ ( i CK )  is  a  monotonical ly  increased 
function  of  ck-  Therefore,  every  point  along  the  imaginary  axis  in 
the  interval  (0,  °°)  corresponds  to  a  single  point  on  the  n-axis  which  lies 
in  the  interval  defined  from  the  intersection  with  the  wave  front  (i.e. 
n  =  vfr  for  the  quasi-longitudinal  and  n  =  /d  for  the  quasi-transverse)  and 
the  origin  (5=  0,  n=  0). 

It  is  interesting  to  examine  what  points  of  the  wave  fields 
correspond  to  the  sides  of  the  branch  cut  from  ©n?  to  infinity.  Direct 
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application  of  Eqs.  (4.16)  yields  the  curves  shown  in  Fig.  6  for  an 
A-I  material  and  Fig.  7  for  a  B-I  material.  The  branch  point  Oq£  is  the 
image  of  point  Aq^  while  the  point  at  infinity  of  the  complex  plane 
corresponds  to  the  origin  of  the  (E»  n)  plane.  These  curves  and  their 
images  in  Fig.  5  are  plotted  with  lines  of  the  same  type  for  easy 
recognition. 

Imagine  now  that  the  quasi-longitudinal  effect  at  a  given  point 
can  be  determined  by  a  sequence  of  values  of  the  parameter  0-j  which 
forms  a  path  on  the  S-j  sheet.  When  the  path  encounters  a  branch  cut, 
the  continuation  of  the  functions  require  a  transfer  from  sheet  to 
$2*  Therefore,  the  curves  traced  on  the  £-n  plane  separate  the  regions 
of  the  wave  fields  which  are  mapped  onto  different  sheets  of  the  Riemann 
surface.  It  should  be  noticed,  however,  that  the  sections  of  the  curves 
located  outside  the  quasi -transverse  wave  front  (see  Fig.  7)  do  not  have 
any  effect  on  the  quasi -transverse  wave  field,  since  that  region  either 
does  not  have  a  representation  or  is  mapped  on  the  real  axis  of  the 
S2  sheet  (see  section  4.2). 

4.3.2  Type  II  Riemann  Surface 

Four  of  the  five  material  cases  satisfy  conditions  of  the  elastic 
parameters  that  are  in  agreement  with  the  inequalities  associated  with 
the  existence  of  the  Type  II  Riemann  surface.  The  study  of  the  Riemann 
surface  will  be  separated  in  two  parts  according  to  whether  the 
quantity  N2  is  positive  or  negative.  This  has  a  direct  effect  on 
whether  the  branch  points  Op  are  located  on  the  $2  or  the  S-|  sheet 
(Appendix  B). 
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4. 3. 2.1  The  Condition  N2  >  0 

The  complete  configuration  of  branch  cuts  in  the  case  under 
examination  is  given  in  Fig.  8.  For  all  real  values  of  6  the  discussion 

K 

in  section  4.3.1  is  directly  applicable.  Every  point  of  the  region  inside 
the  wave  front  is  represented,  in  general ,  by  a  complex  value  of  e  . 
Particular  attention  will  be  given  to  purely  imaginary  values 
specifically  the  part  of  the  imaginary  axis  (0,  ^g2)  which  lies  between 
the  origin  and  the  branch  point  0g2  =  Ug2  on  each  of  the  sheets  and  S2- 
In  that  range  of  values  of  s  the  polynominal  Q  (Uk)  is  a  real  positive 
quantity,  smaller  in  magnitude  than  (b+d+L^)  .  It  becomes  zero  at 
=  ?02‘  num^er  extrema,  determined  from  the  expression 


dfl ( i 4  )  2, 

-JT- s  K^K1K2^  =  0 


(4.18) 


depends  on  whether  M  is  positive  or  negative. 

Let  M<0.  Then  Q  (i^)  has  one  local  maximum  at  ^=0  and  is 
monotonically  decreasing.  The  calculation  of  the  derivative  of  A  (ic  ) 

K  K 

with  respect  to  ;  gives  the  expression 

K 


d  eAu»> 

d^  2bdx ( i c  7 

<  < 


(4.19) 


where 


M?  )  =  L  +  (-D  -1-— 

/xkt 


(4.20) 
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Since  the  Type  II  complex  plane  satisfies  the  conditions  K2  <  0,  N-|  <0 
(Table  B.l),  it  can  be  immediately  seen  that  4»-|(c-|)  >  0,  i.e.  Xi  ( i )  is  a 
monotonically  increasing  function.  Moreover, 


2bN~ 

=  TwT 

(4.21) 

*2^0Z)  =  "  “ 


For  N3  <  0  (i.e.  material  case  D),  ^(i^)  is  negative  in  (0,  Cq?)-  As 

a  result,  X2(is2)  1S  monotonicai ly  decreasing.  The  sides  of  the  branch 

cuts  on  the  imaginary  axis  are  the  images  of  the  closed  curves  through 

the  origin  (see  Fig.  9).  Each  side  of  the  branch  cut  and  the  corresponding 

curves  in  Fig.  9  are  plotted  with  the  same  type  of  line. 

The  foregoing  considerations  lead  to  the  following  interpretations. 

Every  point  of  the  interval  (0,  of  the  imaginary  axis  of  sheet  S-j 

corresponds  to  a  single  point  of  the  interval  (/b^  AQ2)  of  the  n-axis 

(Figs.  8  and  9).  The  points  of  the  quasi-longitudinal  wave  field  between 

n=AQ2  and  n=*F  have  a  representation  on  the  interval  (cQ2,  °)  of  the 

imaginary  axis  of  sheet  S2.  It  was  mentioned  in  section  4.3  that  there 

are  no  values  on  the  complex  plane  to  represent  a  quasi-longitudinal  wave 

field  contained  in  the  triangularly  shaped  area  of  the  quasi -transverse 
★ 

wave  front.  Let  n  denote  the  n-coordinate  of  the  double  point  of  the 

★ 

quasi -transverse  wave  front  and  6  the  corresponding  value  of  parameter  e2 
(see  Table  1).  Then,  every  point  of  the  interval  (/F,  n  )  of  the  quasi- 
transverse  wave  field  has  a  double  representation  on  the  real  axis  of  the 
Sp  sheet  in  the  intervals  (0,  6  )  and  (0,  -6*).  It  can  be  seen  from 
Eqs.  (4.15)  and  (4.16)  that  every  point  of  the  n-axis  in  the  interval 
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1 

1 


★ 

(n  ,  0)  requires  A^(e^)  to  be  a  real  quantity.  The  complex  values  of 
for  which  *K(9K)  becomes  real  form  a  curve  on  the  sheet  S£.  In  fact 

there  are  two  such  curves,  denoted  as  A.  and  BQ  in  Fig.  8,  which  emanate 

0  0 

★  * 

respectively  from  points  -0  and  0  .  The  continuity  of  the  wave  fields 
makes  the  curve  A.  represent  points  of  the  quasi-longitudinal  wave  which 
lie  on  the  interval  (n*»  0)  when  this  interval  is  considered  as  the  limit 
of  the  £  >  0  half  plane.  For  the  quasi-transverse  wave  the  same  limit 

yields  the  curve  B  .  In  effect,  curves  A  and  B  represent  the  boundaries 

0  0  0 

in  the  complex  domain  which  separate  the  regions  of  mapping  the  quasi¬ 
longitudinal  and  quasi-transverse  wave  fields. 

When  condition  Ng  >  0  is  satisfied  (i.e.  material  case  A)  the 
function  2 )  (Eq.  4.20)  changes  sign  from  positive  to  negative  at  a 
point  located  inside  the  interval  (0,  ?q2).  This  point  can  be  determined 
as  the  root  of  the  equation  ^2(c2)  =  Ttle  eliminati°n  of  the  radical 
transforms  this  equation  to 


K  K2c24-2M^  +  l2(b'd-^--  =  0 
w  L  4abd^ 


(4.22) 


The  product  of  the  squares  of  the  two  roots  in  ^  1S  written  as 

[L(b-d)  +  M]  [L(b-d)  -  M]  (4.23) 

4abd2  K] K„ 


and  is  a  negative  quantity  under  the  existing  values  of  M,  l<2  and  Ng.  The 
positive  real  root  in  ^  is 

C*  =  [-  [/id  M+L  /c2^]  /  /ad  K^]  55  (4.24) 


which  can  be  shown  to  satisfy  c2  <  Sq2-  Consequently,  A2(ic2)  is 

★ 

monotonically  increasing  in  the  interval  (0,  c2)  anc'  monotonically 

"At  "ft  "ft 

decreasing  in  (c2,  Cq2).  ^  02  =  ic2*  anc*  n  is  P°"'nt  which  is  mapped 

*  * 
by  e2  in  the  complex  domain,  the  interval  (n  ,  0)  of  the  n-axis  has  a 

double  representation  shown  by  curves  C.  and  D_,  depending  on  whether 

it  is  considered  as  the  limit  of  positive  or  negative  values  of  £.  The 

curves  CQ  and  Da  on  the  Riemann  surface  (Fig.  8)  separate  the  regions  of 

mapping  the  two  types  of  wave  fields  (Fig.  10). 

The  case  M  >  0  is  now  examined.  In  this  case  a  second  local 

extremum  of  the  function  Q(ii^)  exists  in  the  interval  (0,  £Q2)  at  point 


Ck0 


( -  M/  K-j  K2 )  ** 


where  e  g  <  Cg2<  However,  the  general  behavior  of  ( i )  is  not  altered 
and  the  separation  of  the  regions  of  the  complex  domain  corresponding  to 
the  quasi-longitudinal  and  quasi -transverse  wave  fields  is  similar  to  the 
case  >  0,  M  <  0  which  was  discussed  in  this  section.  The  inequalities 
M  >  0  and  <  0  cannot  be  simultaneously  satisfied. 


4. 3. 2. 2  The  Condition  N2  <  0 

The  branching  of  the  Riemann  surface  in  this  case  is  given  in  Fig.  11. 
The  branch  cut  (-1  /  /aT  1  /  /a)  of  the  sheet  represents  the  quasi - 
longitudinal  wave  front  while  the  branch  cuts  (-0qi>  0q^)  of  the  S2  sheet 
and  ( Oqi  ,  1  /  /cf),  (-0Q1,  -1  /  */d)  of  the  sheet  represent  the  quasi- 
transverse  wave  front.  Equations  (4.15)  and  (4.16)  are  employed  again  to 
show  that  real  values  of  in  the  intervals  (1  /  /a”,  1  /  •/cf)  and  (0qi»°°  ) 
represent  points  of  the  quasi-longitudinal  wave  field  in  the  sections 


(/a",  /d)  and  (Cq-j  ,  0)  of  the  £-axis  (Fig.  12).  There  is  no  representation 
of  this  wave  field  for  points  of  the  £-axis  between  £=  /cf  and  £  =  £q-|. 
Correspondingly,  real  values  of  e2  in  the  interval  (oQ1,  00 )  of  The  S2 
sheet  are  the  images  of  points  (£q-j,  0)  of  the  quasi-transverse  wave  field. 
On  the  Riemann  surface,  the  interval  (£q-j  »  >^f)  of  the  £-axis  has  a  double 
representation  which  includes  the  upper  side  of  the  branch  cut  (e^,  1  /  */d) 
on  the  sheet  and  the  section  (0^,  0^)  of  sheet  S2>  Here,  0^  is  the 
point  on  the  S2  plane  corresponding  to  that  half  tangent  plane  which  passes 
through  (£  =  /d,  n  =  0)  and  which  is  not  parallel  to  the  n-axis. 

For  purely  imaginary  values  on  sheets  S1  and  S2  the  discussion 
is  analogous  to  that  given  for  the  case  >  0.  The  separation  of  the 
complex  domain  for  either  >  0  or  <  0  is  given  in  Fig.  11  while  the 
corresponding  wave  fields  for  material  cases  C-II  and  E-II  are  shown  in 
Figs.  12  and  13. 

4.3.3  Type  III  Riemann  Surface 

Figure  14  shows  the  two  sheets  of  the  Type  III  Riemann  surface. 

For  real  values  of  9  the  discussion  in  section  4.3.1  is  directly 

K 

applicable.  Let  eQl  =  and  0q2  =  kq2  be  the  branch  points  on  the 

imaginary  e^-axis.  They  represent  the  images  of  points  AQl  and  AQ2 

of  the  wave  fields  (Figs.  15,  16).  Each  side  of  the  branch  cut 

(Oqi  ,  0q2 )  corresponds  to  a  curve  connecting  points  AQ-]  and  Aq2>  The  side 

of  the  branch  cut  drawn  with  a  dashed  line  (Fig.  14)  corresponds  to  the 

curve  which  is  plotted  with  the  same  type  of  dashed  line  between  points 

Aqi  and  AQ2  (Figs.  15,  16).  It  can  be  shown  that  in  the  intervals 

(0,  cQ2)  and  Uq-|,  »)  of  sheet  S1  the  function  i/»-|(c-|)  is  positive,  thus 
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making  A^(i^^)  a  real  monotonically  increasing  function.  In  the  same 
intervals  of  sheet  S2,  when  N3  >  0,  the  function  Possesses  two  rea^ 

roots  as  a  solution  to  Eq.  (4.22): 


.1*  _  h-L  ^V3^  T-M+L  /-c2N,/ad~ 

‘  ¥5  !  •  and  5  =  —or-1-: 


(4.25) 


i*  2*  1*  1*  2*  2* 

where  £  <  £ .  The  points  e  =  it  and  0  =  it  are  the  images  of 

1*  2* 

points  n  and  n  of  the  n-axis  which  define  a  section  with  a  double 
representation  on  the  complex  domain.  The  curves  A.  and  BD  (Fig.  14) 

o  u 

/  1*  2*x 

map  the  section  (n  ,  n  )  of  the  quasi-longitudinal  wave  field  when  this 
section  is  taken  as  the  limit  of  positive  £  values  (curve  A0)  or 
negative  £  values  (curve  B  ).  For  the  quasi -transverse  wave  field  curve 

C7 

/l*  2*x 

A  maps  the  section  (n  ,  n  )  when  it  is  taken  as  the, limit  of  negative 

0  ‘  ‘  , 

£  values  and  curve  B.  maps  the  same  section  as  the  limit  of  positive  £ 
values.  The  region  enclosed  by  Aa  and  B.  on  the  S0  sheet  (Fig.  14)  maps 

0  0  C 

part  of  the  quasi-longitudinal  wave  field. 

It  is  not  difficult  to  verify  the  following  correspondence  between 


sections  of  the  n-axis  of  the  wave  fields  and  intervals  on  the  Riemann 
surface  (see  Figs.  14  and  15). 


69 


MM 

I 

I 

I 


**?"■ 


mmmrnrnm 


Wave  Field 

Interval  on  n-axis 

Image  on  Riemann  Surface 

Sheet 

(/"b,  A02) 

O 

vT< 

o 

ro 

S1 

(A02,  n1  ) 

l* 

(?02>  ^  ) 

S2 

Quasi-Longitudinal 

(  1*  2*. 

In  »  n  ) 

Curve  A„ 

0 

S2 

(n2  .  Aq.|  ) 

(t2*,  t0)! 

S2 

<Aor  °> 

(?01» 

S1 

</i.  n1*) 

(0,  C1*) 

S2 

Quasi -Transverse 

<n’\  n2*) 

Curve  B 

0 

S2 

(n2*,  0) 

U2  >  °° ) 

S2 

Tho  case  N.,  <  o  (i.e.  material  case  D-III)  provides  one  real  root  of 
Eq.  (4.22)  which  is  larger  in  absolute  value  than  c_, .  Curves  C  and  Da 
of  Fig.  14  are  now  the  boundaries  of  the  regions  of  mapping  the  two  types 
of  wave  fields  which  are  shown  in  Fig.  16. 

4.3.4  Type  IV  Riemann  Surface 

The  branch  points  +0^  and  +0q2  in  this  case  are  located  on  the  real 
8k  axis  (Fig.  17).  When  condition  N2  >  0  is  satisfied  it  can  be  shown 
that  the  quasi-longitudinal  wave  field  is  mapped  entirely  on  the  sheet 
and  the  quasi-transverse  entirely  on  the  S2  sheet  (see  Fig.  18).  The 
same  thing  holds  true  when  <  0  except  for  point  located  on  and  inside 


the  triangularly  shaped  area  of  the  quasi-transverse  wave  adjacent  to 
the  c-axis  which  is  mapped  on  the  branch  cut  (0^  ,  1  /  /cf)  of  the  sheet 
(see  Figs.  19  and  20). 
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5.  NUMERICAL  METHODS  AND  RESULTS 


5. 1  General  Remarks 

The  expressions  for  the  displacements  and  stresses  given  in 
Chapters  2  and  3  determine  the  dynamic  response  in  an  anisotropic  medium 
as  functions  of  the  complex  parameter  e  .  As  time  increases,  the  value 
of  for  a  given  point  traces  a  path  on  the  Riemann  surface.  A 
particular  value  of  on  this  path  determines  the  dynamic  solution  at  the 
corresponding  time.  The  stresses  are  evaluated  by  direct  substitution 
of  in  the  given  expressions,  while  the  displacements  involve  a 
quadrature  in  e^.  The  result  of  these  operations  is,  in  general,  a  complex 
quantity  the  real  part  of  which  represents  the  answer  sought. 

The  principal  objectives  of  this  chapter  are: 

a)  to  focus  attention  on  certain  numerical  problems  which 
appear  in  the  solution  of  anisotropic  wave  propagation 
problems  by  the  method  of  self-similar  potentials  and  to 
explain  the  techniques  which  have  been  used  to  overcome 
them;  and 

b)  to  provide  numerical  results  in  the  form  of  time  histories 
of  displacements  and  stresses  for  the  specific  dynamic 
problems  solved  in  this  study  and  to  discuss  the  most 
characteristic  variations  arising  from  different  classes 
of  anisotropy  as  compared  to  the  isotropic  case. 

Two-dimensional  problems  are  treated  first  in  the  various  anisotropic 
media.  The  elastic  constants  and  density  that  were  used  to  obtain 
numerical  results  were  taken  from  references  [2,  17].  The  paths  of  on 
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sheets  and  S£  of  the  Riemann  surface  are  shown  for  typical  points,  in 
particular  those  located  inside  the  sector  formed  from  the  cusps  and  the 
origin.  Subsequently,  the  numerical  techniques  that  are  applicable  in 
obtaining  the  displacements  for  a  sudden  application  of  a  point  force  in 
transversely  isotropic  materials  are  discussed  to  the  extent  that  they 
differ  from  the  two-dimensional  case. 

The  contributions  of  the  quasi-longitudinal  and  the  quasi-transverse 
wave  fields  at  points  not  near  the  wave  fronts  are  quantities  of  opposite 
sign.  Their  superposition  often  turns  out  to  involve  small  differences 
of  large  quantities.  Therefore,  in  order  to  avoid  erroneous  final 
results,  the  quantities  corresponding  to  individual  waves  must  often  be 
calculated  much  more  accurately  than  might  at  first  be  thought  necessary. 

5.2  Two-Dimensional  Problems 

The  discussion  of  the  two-dimensional  case  concentrates  solely  on 
the  plain  strain  problems  in  orthotropic  materials  excluding  the  much 
simpler  antiplane  problem  which  is  briefly  mentioned  in  section  5.3.2. 

The  methodology  presented  in  the  next  three  sections  not  only  applies  to 
all  the  two-dimensional  problems  of  this  study,  but  also  forms  the  core 
of  the  three-dimensional  case. 

5.2.1  The  Determination  of  eK 

The  determination  of  the  complex  parameter  0^  as  a  function  of  the 
point  coordinates  (x,  y)  and  the  time  variable  t,  given  implicitly  by 
Eq.  (2.16),  requires  the  solution  of  a  complete  fourth  degree  equation 
(see  Eq.  (4.14)).  It  is  convenient  to  use  a  numerical  procedure  for  this 
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calculation. 

As  time  increases,  the  given  point  will  first  have  a  representation 
on  the  Riemann  surface  when  the  tangent  parallel  to  the  horizontal  axis 
passes  through  the  point.  The  corresponding  value  of  is  zero.  Consider 
now  that  at  time  t  the  value  of  the  parameter  is  known  and  given  by 
e".  A  numerical  solution  of  the  differential  equation 

36<  1  1 

IT  =  "  TIFT  =  x  +  y  WFT 

yields  the  new  value  en+1  at  time  t  =  t  +At.  A  predictor- cor rector 
iterative  method  is  used  to  solve  this  equation.  The  number  of  iterations 
required  to  achieve  a  certain  number  of  identical  significant  digits  in 
the  predicted  and  the  corrected  value  of  is  a  criterion  for  deciding 
what  value  of  time  interval  At  should  be  used.  For  example,  when  con¬ 
vergence  is  not  attained  after  a  preset  maximum  number  of  iterations,  it 
is  an  indication  that  the  rate  of  change  of  is  very  large  and  it  is 
convenient  to  reduce  the  time  interval.  It  is  reduced  by  a  factor  of  2 
in  the  program.  On  the  other  hand,  if  the  convergence  criterion  is 
satisfied  in  less  than  a  preset  minimum  number  of  iterations,  the  rate  of 
change  of  is  very  small  and  the  time  interval  is  increased  by  a  factor 
of  2.  In  this  process,  in  order  to  avoid  excessively  large  time  intervals, 
an  upper  limit  of  At  is  imposed.  However,  if  e  is  approaching  a  singular 

K 

point.  At  would  assume  values  smaller  than  any  preset  lower  limit.  In 
the  neighborhood  of  such  a  singularity  the  values  of  are  obtained  by  a 
modified  method  to  be  explained  in  the  next  section. 


i 


n  1  1 

Since  the  evaluation  of  e"  depends  on  the  previously  calculated 
value  e",  the  errors  of  approximation  accumulate  causing  an  increasing 
drift  of  the  calculated  values  from  the  exact  ones.  To  eliminate  this 
effect,  the  use  of  the  Newton-Raphson  method  [36]  with  a  higher 
convergence  criterion  is  employed  at  selected  time  instances.  Thus,  the 
desired  values  are  evaluated  with  better  accuracy  and  the  accumulation 
of  error  is  canceled. 

The  combined  use  of  the  predictor-corrector  and  the  Newton-Raphson 
iterative  methods  provides  a  very  efficient  technique  for  the  numerical 
solution  of  Eq.  (2.16).  Direct  application  of  only  the  latter  method 
could  fail  to  converge  or  even  might  converge  to  the  wrong  value  if  the 
initial  guess  was  not  good  enough.  The  predictor-corrector  method 
provides  a  fast  means  of  obtaining  a  sufficiently  good  approximation  in 
addition  to  the  criterion  for  an  approaching  singularity. 

Particular  care  should  be  exercised  in  the  computation  of  the 
numerical  value  of  the  function  x^(e^).  The  double  radical  expression  of 
Eq.  (2.19)  is  calculated  in  two  steps.  In  the  first  step,  the  value  of 
the  square  root  of  the  polynomial  Q(e  )  is  computed  by  means  of  a  special 

1C 

routine  as  shown  in  Appendix  D.l.  Since  the  term  (-1)K  assumes  a  definite 
value  on  each  of  sheets  and  S2  of  the  Riemann  surface,  every  value  of 
yields  a  unique  value  for  the  expression  under  the  outside  radical  on 
the  right  hand  side  of  Eq.  (2.19).  In  the  second  step,  the  calculation  of 
the  square  root  of  this  last  quantity  is  now  straightforward  (see  Appendix 


5.2.2  The  Behavior  Near  the  Wave  Fronts 

In  Chapter  4  it  has  been  mentioned  that  the  condition  satisfied  by 
points  on  the  wave  fronts  can  be  expressed  by  the  equation 

«;<%>  *  -x  -y  -  °  <5-2> 


where  e^0  is  the  value  of  the  parameter  corresponding  to  point  (x,  y) 
on  the  wave  front  at  time  tg.  The  differential  equation  (5.1)  becomes 
singular  at  this  point.  At  times  very  close  to  tQ  the  rate  of  change  of 
0^  is  very  large.  As  a  result,  the  time  interval  At  required  for  con¬ 
vergence  in  the  predictor-corrector  method  becomes  very  small.  In  such  a 
case  it  is  convenient  to  use  a  modified  version  of  the  method  to  obtain  a 
better  approximation  of  the  predicted  value  by  taking  account  of  the 
behavior  of  Eq.  (5.1)  at  points  near  the  wave  fronts. 

Consider  the  power  series  expansion  of  the  function  6^(0^)  about 

the  value  e  n,  i.e. : 
kU 


6"1  (0  ) 

W '  6k<%>  *  fV'UK'Vo* + 


(5.3) 


Since  fi^(e  q)  =  0  it  can  be  rewritten  as  follows: 


6' (0  )  =  (0 
K  K  K 


[-yA^ 


(<w-  'vw y 


2! 


(5.4) 


Equation  (5.1)  can  now  be  integrated  between  tQ  and  t  to  give 
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This  last  equation  suggests  that  for  points  near  the  wave  fronts  the 
change  of  the  parameter  ©  can  be  closely  approximated  as  being 
proportional  to  the  square  root  of  the  change  of  time.  Based  on  this 
result,  the  right  part  of  Eq.  (5.1)  can  be  expressed  approximately  as  a 
linear  function  in  time  divided  by  the  term  /t-tg,  which  gives  rise  to  the 
singularity  at  t=tg.  The  evaluation  of  the  integral  of  this  quantity  is  a 
simple  matter  (see  Appendix  D.3)  and  provides  a  very  good  first  estimate 
of  e  .  There  is  need  for  only  a  small  number  of  subsequent  iterations  u> 
determine  e  to  the  desired  level  of  accuracy. 

K 

The  application  of  the  Newton- Raphson  method  at  points  near  the 
wave  fronts  is  also  ineffective.  The  reason  is  the  existence  of  a  double 


root  of  the  equation  =  0  at  e^Q,  which  makes  it  almost  impossible  to 

6  (en)  * 

apply  the  correction  term  — 77“rrr  to  go  from  iteration  to  iteration.  The 

6  (0  ) 

K  '  K 


problem  can  be  bypassed  by  seeking  the  root  of  a  new  function 
6  ( ®  ) 

g(eK>  =  ^V0K)  since  e  Q  is  a  single  root  of  g(e  )  [36,  section  8.6]. 


5.2.3  The  Description  of  the  Paths  of  0^ 

Let  the  coordinates  (x^,  y^)  define  a  point  on  the  x>0,  y>0  quarter 
plane  for  which  displacements  and  stresses  are  to  be  found.  From  the 
discussion  of  Chapter  4  it  becomes  obvious  that,  for  a  given  material,  the 
path  of  the  complex  variable  0  traced  on  the  complex  domain  as  time 
increases  depends  entirely  on  the  ratio  Y=-^.  Let  t„  denote  the  arrival 
time  of  the  quasi-longitudinal  wave  and  t$  the  arrival  time  of  the  inside 
convex  part  of  the  quasi-transverse  wave  front.  For  points  inside  the 
sector  formed  from  the  cusps  and  the  origin,  t$-|  indicates  the  arrival  time 
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of  the  external,  concave  part  of  the  quasi -transverse  wave  front  between 
the  two  cusps  and  t$2  the  arrival  time  of  the  section  between  one  of  the 
cusps  and  the  double  point  B^E  of  the  same  front  (Fig.  4).  Consider 
now  the  tangents,  parallel  to  the  x-axis  at  the  intersection  of  each  wave 
front  with  the  positive  axis.  Let  tpQ  =  yR  /  /F and  t$0  =  yR  /  /d~ be  the 
corresponding  times  that  these  tangents  pass  through  point  (xR,  yR).  It 
is  always  true  that  tpQ  <_  tp<ts]i  tS2-ts  and  tsO-ts’ 

Since  different  values  of  Y  correspond  to  points  with  different 
locations  relative  to  the  cusps  (Fig.  7),  the  shape  of  the  curves  that 
represent  the  paths  of  e  varies  accordingly.  A  number  of  such  curves  is 
shown  in  Figs.  21-27  for  different  material  cases  and  different  values  of 
Y.  It  is  of  interest  to  observe  that  for  values  of  y  corresponding  to 
points  inside  the  sector  included  between  the  cusps  and  the  origin  and 
for  times  after  the  arrival  of  the  first  wave,  the  complex  value  of 
which  represents  the  quasi -longitudinal  wave  field  eventually  approaches 
a  real  value  as  time  approaches  tsl  (Figs.  22,  25,  27).  Moreover,  at  times 
after  ts2  the  parameter  assumes  complex  values  again.  It  is  not 
difficult  to  show  that  at  these  two  instants  Eq.  (5.1)  is  singular  while 
at  times  just  before  t$^  or  just  after  tg2  the  behavior  is  as  'scribed 
in  section  5.2.2.  The  conclusion  that  seems  to  suggest  itself  is  that 
the  boundary  of  the  triangularly  shaped  area  of  the  cusped  front  is  not 
only  part  of  the  quasi-transverse  wave  front  but  is  also  a  front  of  the 
quasi-longitudinal  wave.  The  enclosed  area  does  not  experience  any 
disturbances  of  the  quasi-longitudinal  type. 


5.2.4  The  Computation  of  Displacements  and  Stresses 

Selected  time  histories  of  displacements  and  stresses  are  given 
for  the  following  two-dimensional  problems  (Fig.  28): 

a)  Horizontal  load  in  a  full-space  ( FS-H ) 

b)  Vertical  load  on  the  surface  of  a  half-space  (HS-V) 

c)  Horizontal  load  in  the  interior  of  a  half-space  (HS-IH) 

The  load  consists  of  a  suddenly  applied  line  force.  It  is  instructive  to 
include  in  certain  cases  the  displacements  for  the  impulsive  application  of 
the  load.  At  every  time  t  the  corresponding  value  of  the  parameter 
e”  on  the  path  described  is  substituted  in  a  straightforward  manner  to 
calculate  the  values  of  the  stresses  for  the  suddenly  applied  force  and 
the  two  displacements  for  the  impulsive  load.  The  displacements  for  the 
suddenly  applied  force  are  found  by  performing  a  quadrature  along  a  path 
with  end  points  e  =0  and  e  =en.  Any  path  between  these  two  points  can  be 
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used  provided  the  conditions  for  the  Cauchy-Goursat  theorem  [9]  are 
fulfilled.  However,  it  is  convenient  to  use  the  path  that  traces  on 
the  complex  domain  as  time  increases.  Then,  the  displacements  at  time  Vl 
are  computed  from  the  values  at  time  t  and  the  contribution  of  the  integral 
between  values  e11  and  en+1.  The  latter  may  be  obtained  by  use  of  any 
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numerical  integration  scheme.  In  the  present  study  satisfactory  answers 
have  been  found  from  the  application  of  the  trapezoidal  rule  of  integration. 
This  may  be  due  to  the  fact  that  the  time  interval,  which  in  effect  sets 
the  integration  interval,  has  been  proportioned  in  accordance  with  the 
rate  of  change  of  (and  all  functions  of  e^)  with  time. 
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All  results  are  plotted  against  the  dimensionless  quantity  t/t$. 

It  can  be  seen  that  in  all  cases  the  displacement  in  the  direction  of 
the  suddenly  applied  force  is  not  bounded.  The  displacement  in  the 
direction  perpendicular  to  the  loading,  as  well  as  the  stresses, 
asymptotically  approach  constant  values  for  large  times.  For  the  vertical 
force  on  the  surface  of  the  half  space  these  constant  values  are  in 
agreement  with  the  static  solutions  of  reference  [20].  Solutions  to 
dynamic  problems  in  isotropic  materials  can  be  obtained  from  the  expressions 
of  Chapter  2  by  imposing  the  conditions  of  Appendix  G.  For  sufficiently 
large  values  of  the  ratio  (t/ts)  the  dynamic  quantities  are,  as  expected, 
in  agreement  with  the  static  solutions  of  references  [21,  22,  24]. 

The  basic  features  of  the  dynamic  response  at  a  point  of  an 
orthotropic  material  can  best  be  exemplified  by  considering  the  full 
space  problem  where  the  effects  from  the  head  and  Rayleigh  waves  are  absent. 
Let  rc  indicate  any  value  of  y  such  that  the  corresponding  point  en¬ 
counters  the  concave  part  of  the  cuspidal  front  and  let  Tg  be  any  other 
value  of  y .  Consider  also  an  anisotropic  medium  of  the  classes  A-II,  B-I 
or  C-II.  The  time  histories  for  points  with  y=0.5,  y=0.8  and  y=10  are 

given  in  Figs.  29  to  33.  Whenever  Y=rg,  the  arrival  of  each  of  the  two 

★ 

wave  fronts  is  clearly  apparent  from  the  resulting  singularities  in  the 
diagrams  of  stresses  as  well  as  in  the  diagrams  of  the  displacements 
caused  by  the  impulsive  load.  When  Y=fc  there  are  four  such  singularities 
corresponding  to  the  four  encounters  with  sections  of  the  wave  fronts. 

While  the  point  under  consideration  is  inside  the  triangularly  shaped 

it 

In  the  diagrams  of  the  time  histories  of  the  displacements  and  stresses, 
these  singularities  are  denoted  by  an  arrow. 
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area  or  "lacuna"  of  the  quasi -transverse  wave  front  no  elastic  vibrations 
are  experienced,  as  can  be  seen  from  the  vanishing  of  the  stresses 
(Figs.  31b  and  33b)  and  the  displacements  for  the  impulsive  load 
(Figs.  31c  and  33c).  The  effect  on  the  displacements  for  the  step  force 
is  best  described  as  "zero  change"  shown  as  the  flat  part  of  the  response 
in  Figs.  31a  and  33a.  The  comparison  of  the  time  histories  for  the  same 
value  of  y  but  for  different  materials  shows  that  any  differences  occur 
for  t/t  <1.  It  is  interesting  to  point  out  the  rapid  change  in  the 
response  before  the  arrival  of  the  quasi-transverse  wave  shown  in 
Figs.  30b  and  30c.  Such  a  change  could  be  confused  with  the  arrival  of 
the  second  wave  front,  but  it  is  merely  caused  by  the  nearness  to  the 
singularity  at  the  cusp.  If  the  response  is  parametrically  examined  in 
terms  of  Y  it  can  be  seen  that  this  change  becomes  more  pronounced  as  Y 
increases  and  eventually  evolves  into  the  two  singularities  that  correspond 
to  the  encounter  of  the  boundary  of  the  lacuna,  which  was  discussed 
previously.  All  classes  of  plane  strain  orthotropic  media  with  cuspidal 
fronts  exhibit  this  kind  of  behavior  in  the  time  histories  of  the  stresses. 

On  the  other  hand,  the  dynamic  response  of  materials  with  strictly  convex 
fronts  (Fig.  29)  is  very  similar  to  the  response  observed  in  isotropic  media. 

The  half  space  solutions  where  the  load  is  applied  on  the  surface  are 
presented  in  Figs.  34-38  for  the  same  classes  of  orthotropic  materials  as 
in  the  previous  problem  (i.e.  A-II,  B-I  and  C-II).  Displacements  and 
stresses  are  given  for  points  with  y=  0.8,  y=  2.2  and  y=  10.0.  It  is  a 
well  known  fact  that  the  presence  of  a  free  boundary  causes  the  generation 
of  head  waves  and  surface  (Rayleigh)  waves  which  may  influence  the  dynamic 
response  depending  on  the  location  of  the  point.  The  head  wave  regions  for 
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materials  with  strictly  convex  fronts  are  defined  in  exactly  the  same  way 
as  in  the  case  of  isotropy.  The  presence  of  cusps  alters  the  geometry  of 
these  regions  as  can  be  shown  by  the  shaded  areas  in  Fig.  39.  It  is 
interesting  to  remark  that  for  materials  in  the  classes  C-II,  E-II  and 
E-IV,  which  have  lacunas  crossed  by  the  x-axis,  there  exists  a  leading  and 
a  trailing  wave  front  (Fig.  39(b)).  This  is  a  result  of  the  absence  of 
any  quasi-longitudinal  type  disturbances  in  the  section  of  the  y  =  0  plane 
contained  inside  the  lacuna  of  the  quasi-transverse  wave  front,  thus  making 
the  existence  of  head  wave  disturbances  unnecessary.  The  effect  of  the 
Rayleigh  waves  becomes  quite  noticeable,  as  expected,  for  points  close  to 
the  free  surface  (i.e.  for  large  values  of  y).  This  effect  is  shown  in 
Fig.  38  by  the  rapid  change  in  the  response  just  after  the  arrival  of  the 
quasi-transverse  wave. 

The  time  histories  of  the  surface  displacements  uxg,  uyQ  and  the 
surface  stress  axQ  for  the  third  problem  of  a  subsurface  horizontal  force 
are  given  in  Figs.  40-44.  In  this  case,  a  head  wave  is  present  in  order 
to  satisfy  boundary  conditions  on  the  free  surface  prior  to  the  arrival  of 
the  reflected  quasi-transverse  and  after  the  passage  of  the  reflected 
quasi-longitudinal  wave.  Here,  both  of  these  two  reflected  waves  are 
generated  after  an  incident  quasi-transverse  wave  encounters  the  y  =  0 
plane.  The  presence  of  cusps  in  the  incident  wave  front  will,  of  course, 
complicate  the  geometry  of  the  head  wave  regions.  However,  any  con¬ 
tribution  that  may  affect  the  response  is  included  automatically  in  the 
general  solutions  by  the  Smirnov-Sobolev  method  and  consequently  in  the 
given  time  histories.  The  effect  of  the  Rayleigh  waves  becomes  more 
pronounced  at  points  sufficiently  far  from  the  origin  (Fig.  44). 
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5.3  Three-Dimensional  Problems 

The  numerical  solution  of  the  problem  of  the  sudden  application  of  a 
point  force  on  the  surface  of  a  transversely  isotropic  medium  incorporates 
the  same  techniques  used  for  the  two-dimensional  plane  strain  problems, 
up  to  the  point  that  the  path  of  0^  is  defined  on  the  Riemann  surface.  In 
order  to  show  how  the  antiplane  problem  is  treated,  the  more  complicated 
case  of  a  tangential  force  will  be  considered  instead  of  the  simpler 
axi symmetric  case  of  a  normal  force. 

5.3.1  The  Contours  of  Integration  CK 

The  formal  solution  of  the  dynamic  Cerruti  problem  for  a  transversely 
isotropic  medium  is  given  by  Eqs.  (3.36)  which  were  derived  from  Eqs. 

(3.26)  by  changing  the  variable  of  integration  from  W  to  e^.  As  W  varies 
in  the  closed  interval  (0,tt)  the  variable  assumes  complex  values  which 
satisfy  Eq.  (3.22).  These  values  trace  a  contour  on  the  Riemann  surface 
along  which  the  integration  indicated  by  Eq.  (3.3F‘  must  be  performed.  The 
end  points  of  the  contour  are  determined  as  the  solutions  of  Eq.  (3.22) 
for  values  W  =  0  and  W  =  ir.  Every  value  of  0  on  the  contour  C  determines 
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the  contribution  to  the  three-dimensional  wave  field  of  a  two-dimensional 
wave  field  located  at  the  corresponding  angular  position  W.  The  numerical 
computation  of  the  values  0^  on  contours  C^  is  performed  in  exactly  the  same 
way  as  in  two-dimensional  problems  since  for  every  value  of  W  an 
equivalent  x  may  be  taken  as  x  =  rcos  W. 

5.3.2  The  Numerical  Treatment  of  the  Three-Dimensional  Solutions 

The  considerations  of  the  previous  section  suggest  that  the  approach 
to  be  followed  in  the  numerical  treatment  of  the  three-dimensional 
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solutions  may  be  divided  into  three  steps. 

In  the  first  step,  the  paths  of  for  the  plane  strain  and  antiplane 
problems  at  the  angular  position  W=0  are  determined  as  time  increases. 

The  numerical  solution  of  Eq.  (3.22)  is  required  to  achieve  this  task  for 
the  plane  strain  problem.  However,  the  complex  parameter  e-j,  which 
defines  the  contribution  of  the  anti  plane  problem,  can  be  explicitly 
expressed  as  a  function  of  x,  y  and  t,  i.e. 

e  tx  ±  y  /!/? /(xi±.gZz.d]/^?  (5.6) 

3  (xZ  +  eyZ/d) 

The  sign  of  the  y-term  is  chosen  to  correspond  with  the  mapping  of  the 

* 

y  _>  0  half  plane  to  the  Ime^  _>  0  complex  plane. 

In  the  second  step,  the  effect  of  t.he  rotational  superposition  on 
the  displacement  field  at  every  value  of  time  is  found  by  performing  the 
quadratures  along  the  contours  (<*1,  2,  3).  It  is  important  to 

determine  whether  the  contours  and  C2  are  composed  of  sections  that  lie 

entirely  on  a  single  sheet  of  the  Riemann  surface  or  not.  The  reason  is 
that  in  the  numerical  computation  of  the  quadratures  it  is  often  convenient 
by  use  of  the  Cauchy-Goursat  theorem  [9]  to  deform  the  actual  contour  into 
a  number  of  straight  line  segments.  Contour  C3  lies  entirely  on  the  02- 
complex  plane.  Notice  also  that  the  singularity  of  the  integrand 
expressions  of  Eqs.  (3.36)  at  =  0  must  be  avoided  by  the  deformed 
contours. 

Finally  in  the  third  step,  a  numerical  integration  scheme  is  used 
to  perform  the  quadratures  over  suitably  deformed  contours  CK-  The 

*The  03-complex  plane  has  a  branch  cut  determined  from  the  function  A3 (03) 
(see  Eq.  (2.32)).  In  this  studyit  is  taken  along  the  real  axis 
connecting  branch  points  +  1  //e. 


displacements  are  calculated  by  adding  the  contributions  of  the  quasi¬ 
longitudinal  wave  (contour  C^),the  quasi-transverse  wave  (contour  C2) 
and  the  shear  wave  of  the  antiplane  problem  (contour  C^).  In  this 
study  Simpson's  rule  is  used  except  for  a  short  segment  at  each  end  of 
the  contour.  There,  the  singularities  introduced  by  the  function  G^e^) 
require  special  consideration,  which  is  the  subject  of  the  next  section. 

5.3.3  The  Behavior  at  the  Ends  of  the  Contours  CK 

Let  the  function  G^e^)  (see  Eq.  (3.25))  be  written  in  the  form 

h  h 

W  =  [reK_t  +  y  tr0K  +  t-y  (5-7) 


The  first  radical  in  the  product  contains  the  left-hand  side  of  Eq.  (3.22) 
for  W  =  0  and  the  second  radical  contains  the  same  expression  for  W  =  u. 
That  is,  the  ends  of  the  contour  C  are  the  zeros  of  the  function  G  ( e  ) . 
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If  indicate  the  initial  point  of  the  contour  it  is  not  difficult  to 
show  that  a  power  series  expansion  about  e*  transforms  Eq.  (5.7)  to  the 


expression 


h 

G  (e  )  =  (e  -e11)  G  (e  ) 

k'k*  'k 


(5.8) 


where  G^e*)  ^  0.  As  a  result,  the  integrands  of  Eqs.  (3.36)  can  be 

0  l< 

written  in  terms  of  the  radical  (e  -e  )  which  behaves  nonanalytically  at 
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and  a  function  which  is  analytic  at  e^.  The  latter  can  be  approximated 
by  its  linear  part  and  integrated  as  shown  in  sections  4.1  and  4.2  of 
Appendix  D  to  obtain  an  integration  formula.  This  formula  is  used  for 
the  numerical  computation  of  the  quadrature  along  a  small  section  at  the 
ends  of  the  contours. 
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If  the  particular  value  of  the  time  variable  corresponds  to  the 
passage  of  one  of  the  wave  fronts  at  the  point  under  consideration,  then 
it  can  be  shown  that  the  function  G  (e  )  can  be  written  as 

K  K 

G  (0  )  =  (0  -e£)  G*(e  )  (5.9) 

K  K  K  <  K 
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where  GK(e^)  t  0.  The  same  procedure  as  previously  is  used  to  develop 
the  integration  formula  given  in  Appendix  D.4.3. 


5.3.4  Sample  Problems 

The  methodology  described  in  section  5.3.2  has  been  used  to 

calculate  the  time  histories  of  the  displacements  u„,  u  ,  u„  when  the 

r  w  y 

applied  load  consists  of  a  suddenly  applied  tangential  point  force  at 
the  origin.  The  wave-front  pattern  generated  for  a  transversely 
isotropic  material  of  class  B-I*  is  shown  in  Fig.  45.  Sample  dynamic 
responses  for  various  values  of  Y**  in  transversely  isotropic  media  of 
classes  B-I  and  A-I  are  given  in  Figs.  46-48. 

The  arrival  of  the  quasi -longitudinal  wave  results  in  a  step 
discontinuity  of  the  radial  and  vertical  displacements  while  the 
circumferential  displacement  behaves  continuously  starting  from  a  value 
of  zero.  A  similar  step  discontinuity  in  the  radial  and  vertical 
displacements  appears  also  upon  the  arrival  of  the  quasi-transverse 
wave  front,  while  the  arrival  of  the  shear  type  (antiplane)  wave  results 
in  a  step  discontinuity  of  the  circumferencial  displacement.  At  large 


^Transversely  isotropic  materials  are  designated  in  this  study  by  the  class 
which  characterizes  the  corresponding  plane-strain  orthotropic  material. 


** 


In  three-dimensional  problems  Y=rp/yR,  where  (rp,  u>r,  yR)  are  the 
cylindrical  coordinates  of  the  point  under  consideration. 
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times  the  displacements  approach  the  static  solutions  found  in  reference 
[48].  When  the  conditions  of  Appendix  G  are  applied,  the  time  histories 
of  the  displacements  correspond  to  the  solutions  for  an  Isotropic  material. 
These  solutions  at  large  times  approach  the  static  solutions  of  reference 
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6.  SUMMARY ,  CONCLUSIONS  AND  RECOMMENDATIONS 
FOR  FURTHER  STUDY 

6.1  Suntnary 

The  method  of  self-similar  potentials  or  functional  invariant 
solutions  was  used  in  this  study  to  obtain  solutions  for  certain  two-  and 
three-dimensional  wave  propagation  problems  in  anisotropic  media.  The 
solutions  were  expressed  in  terms  of  analytic  functions  which  were 
determined  from  the  given  boundary  conditions  in  a  straightforward  manner. 

In  the  process  only  simple  algebraic  operations  were  required,  a  distinct 
advantage  over  alternative  methods  employing  the  more  complicated 
techniques  of  multiple  inverse  transforms. 

The  two-dimensional  cases  solved  wave  propagation  problems  for  the 
sudden  application  of  a  line  force  in  an  orthotropic  medium  under  plane 
strain  conditions.  Each  of  the  propagating  waves  was  parametrized  by  a 
complex  variable  which  was  used  as  the  argument  of  a  corresponding  potential 
function.  Thus,  the  problem  of  reflection  was  easily  handled  by  forming  a 
system  of  linear  simultaneous  equations,  equal  in  number  to  the  boundary 
conditions,  and  solving  for  the  unknown  potentials.  Any  singularities  in 
the  wave  fields  could  be  inferred  from  the  examination  of  the  singularities 
of  the  associated  potential  functions. 

Solutions  of  appropriately  defined  two-dimensional  problems  were 
superimposed  in  a  rotational  manner  to  obtain  the  wave  field  in  transversely 
isotropic  media.  The  displacements  for  the  axi symmetric  problem  of  a  normal 
force  and  the  non-axi symmetric  problem  of  a  tangential  force  on  the  y  =  0 
surface  were  given  in  the  form  of  simple  quadratures.  The  technique  of 
rotational  superposition  made  it  possible  to  solve  these  three-dimensional 
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problems  with  little  additional  difficulty  than  the  two-dimensional  cases. 

The  correct  representation  of  the  wave  fields  in  the  complex  domain 
proved  to  be  extremely  important  in  the  numerical  treatment  of  the 
solutions  and,  therefore,  was  thoroughly  examined.  It  was  shown  that  wave 
fields  in  plane  strain  othotropic  materials  can  be  mapped  onto  four 
different  types  of  Riemann  surfaces.  In  the  numerical  calculation  of  the 
displacements  and  stresses  efficient  techniques  were  used  to  locate  the 
existing  singularities.  The  quadratures  of  the  displacements  in 
transversely  isotropic  media  were  performed  by  using  Simpson's  rule  except 
for  small  segments  at  the  end  points  of  the  integration  contours.  There, 
the  dominant  nonanalytic  behavior  was  taken  into  account  by  employing 
special  integration  formulas. 

6.2  Conclusions 

The  development  of  the  solutions,  the  discussion  of  the  mapping  and 
the  presentation  of  the  results  lead  to  the  following  conclusions: 

1)  The  method  of  the  self-similar  potentials  can  also  be  used 
to  provide  a  straightforward  approach  in  the  solution  of  wave 
propagation  problems  in  certain  anisotropic  materials. 

Isotropy  can  be  treated  as  a  special  case;  however,  the 
computational  effort  to  develop  the  more  general  anisotropic 
solution  is  significantly  larger.  Of  particular  importance 
is  the  fact  that  the  quasi-longitudinal  and  quasi-transverse 
wave  fields  are  mapped,  in  general,  in  regions  which  do  not 
lie  on  a  single  sheet  of  the  Riemann  surface. 
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2)  Under  certain  conditions  on  the  elastic  parameters  the  quasi- 
transverse  wave  front  possesses  cusps  and  double  points  which 
form  triangularly  shaped  lacunas.  Elastic  stresses  do  not 
exist  for  points  inside  the  lacunas  except  when  the  lacuna  is  in 
a  head  wave  region.  In  a  sense,  then,  the  sections  of  the 
quasi- transverse  wave  front  which  form  the  boundary  of  the 
lacunas  can  also  be  thought  of  as  being  part  of  a  quasi¬ 
longitudinal  wave  front  even  though  they  are  moving  out  at 
quasi-transverse  speeds.  As  a  result,  points  which  are  located 
so  that  they  experience  four  encounters  with  various  sections 

of  the  wave  fronts  present  four  singularities  in  the  stress 
field.  At  all  other  point  locations  only  two  such  encounters 
occur.  However,  rapid  changes  in  the  stress  responses  are 
observed  which  are  affected  to  a  varying  degree  by  the  nearness 
to  any  of  the  cusps.  The  behavior  of  materials  with  no  cusps 
in  the  wave  fronts  is  very  similar  to  the  isotropic  case. 

3)  The  solutions  of  the  half-space  problems  demonstrate  the 
existence  of  head  wave  regions  which,  in  the  case  of  materials 
with  cusps  in  the  fronts,  are  somewhat  different  than  those  in 
isotropic  materials.  In  certain  cases  there  exists  a  trailing 
as  well  as  a  leading  head  wave  front  in  contrast  to  a  single 
front  which  is  formed  when  the  quasi-transverse  wave  does  not 
have  cusps. 

4)  The  extension  of  the  method  of  rotational  superposition  permits 
wave  propagation  problems  in  transversely  isotropic  media  to  be 
treated  almost  as  easily  as  any  of  the  two-dimensional  problems. 


It  should  be  emphasized  that  in  this  case  solutions  for  two 
different  materials  are  superposed  to  obtain  the  displacement 
field  in  a  transversely  isotropic  medium. 

6.3  Reconmendations  for  Further  Study 

The  solutions  which  have  been  developed  in  this  study  can  be  used 
directly  for  a  more  detailed  examination  of  the  dynamic  behavior  of 
anisotropic  media.  For  example,  a  parametric  study  of  the  time  histories 
where  the  elastic  constants  vary  continuously  may  make  possible  the 
recognition  of  the  characteristics,  if  any  exist,  in  the  response  that 
may  be  associated  with  the  conditions  for  the  initial  formation  of  the 
cusps. 

A  more  general  time  variation  of  the  applied  force  in  the  two- 
dimensional  problems  can  be  constructed  by  use  of  a  convolution  integral. 
In  such  a  process  solutions  for  different  magnitudes  of  an  impulsive  load 
are  superposed  at  different  times.  In  seismology  and  in  fracture 
mechanics  the  problem  of  an  expanding  dislocation  in  the  interior  of  a 
half-space  is  of  great  importance  and  can  be  solved  for  anisotropic 
materials  by  the  methodology  of  Chapter  2. 

Another  interesting  application  concerns  the  development  of 
solutions  for  two  bonded  half-spaces  each  representing  a  different 
anisotropic  material.  The  generated  reflected  and  refracted  waves  are 
associated  with  potentials  which  can  be  easily  determined  from  the 
boundary  conditions.  The  solution  to  this  problem  may  be  the  first  step 
in  formulating  a  more  general  layered  anisotropic  space. 
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A  challenging  question  is  posed  as  to  whether  a  modified  Smirnov- 
Sobolev  method  could  be  developed  to  solve  the  problems  worked  out  in 
this  study  for  a  more  general  anisotropic  material. 

In  transversely  isotropic  media  the  problems  of  a  point  force,  a 
double  force  with  moment  or  a  point  dislocation  in  the  interior  of  a 
half  space  can  be  obtained  from  the  solutions  of  Chapters  2  and  3  with 
little  additional  effort.  The  case  of  two  bonded  transversely  isotropic 
half-spaces  can  also  be  treated  by  the  technique  of  rotational  super¬ 
position.  However,  the  latter  may  require  some  modification  before  it 
can  be  applied  to  more  general  types  of  anisotropy  than  the  hexagonal 
symmetry.  This  point  still  needs  further  exploration. 
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TABLE  I 


RELATIONS  BETWEEN  ELASTIC  CONSTANTS  AND  THE  SHAPE 
OF  THE  QUASI-TRANSVERSE  WAVE  FRONT 


Note  I 


Note  2. 


1  2*  ^ 

e  ,  e  and  0  are  roots  of  function  ( 0 )  in  (0,  1/  *  d) 
o'**  is  the  root  of  02(e)  in  (0,  g-j ) 

”0  is  the  root  of  D1 (e)  in  (Og1 ,  1/  /d) 

,2*  f  M-L  /-c^N-j/ad  1 


** 

0  =  0 


K1K2 


1/2 


Arrow  indicates  ronotonically  increasing  (/)  or  monotonical ly 
decreasing  (H)  values  in  the  specified  interval. 
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FIG.  3  THE  REPRESENTATION  OF  THE  QUASI-LONGITUDINAL  WAVE  ON  SHEET  S] 
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4  THE  REPRESENTATION  OF  A  CUSPED  QUASI-TRANSVERSE  WAVE  ON  SHEET 
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FIG.  6  WAVE  FIELDS  FOR  MATERIAL  CASE  A- I 
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FIG.  7  WAVE  FIELDS  FOR  MATERIAL  CASE  B-I 
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FIG.  10  WAVE  FIELDS  FOR  MATERIAL  CASE  A-II 


U 

0 


Quasi- Longitudinal 
Wave  Front  — ^ 


Quasi-Longitudinal 
Wave  Front 


112 


FIG.  16  WAVE  FIELDS  FOR  MATERIAL  CASE  D-III 


FIG.  18  WAVE  FIELDS  FOR  MATERIAL  CASE  A-IV 
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FIG.  23  PATHS  OF  FOR  CASE  B-I  AND  y=2.2 
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FIG.  25  PATHS  OF  e  FOR  CASE  C-II  AND  y=0.1 
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FIG.  31c  DISPLACEMENTS  FOR  FS-H  IN  A  CLASS  B-I  MATERIAL  WITH  Y«  0.8 
AND  IMPULSIVE  APPLICATION  OF  THE  LOAD 
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IISPLACEMENTS  FOR  HS-V 


IN  A  CLASS  B-I  MATERIAL  WITH  >=0.8 


149 


0  » 


a)  MATERIAL  CASE  B-I 


b)  MATERIAL  CASE  C-It 


FIG.  39  THE  HEAD  WAVE  REGIONS 
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FIG.  48  DISPLACEMENTS  FOR  A  SURFACE  HORIZONTAL  FORCE  IN  A  CLASS  B-I 
TRANSVERSELY  ISOTROPIC  MATERIAL  WITH  V  =  0.8 


APPENDIX  A 


THE  DERIVATIVES  OF  9 

K 

The  first  and  second  derivatives  of  0^  with  respect  to  x,  y  and  t 
can  be  obtained  from  Eq.  (2.16)  by  following  the  rules  of  differentiation 
for  implicit  functions  [42].  They  are  given  by  the  following  expressions 
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where  6^  is  given  by  Eq.  (2.37). 

Let  f(0  )  and  f'(0  )  denote  any  analytic  function  of  e  and  its 

K  K  K 

derivative  with  respect  to  the  same  variable.  The  first,  second  and  third 
derivatives  of  f(0K)  with  respect  to  x,  y  and  t  are 
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where  p  +  q  +  r  =  3 

In  the  case  of  the  reflected  waves  the  same  formulas  can  be  used.  The 
only  differences  are  that  x  (e  )  must  be  replaced  by  -X  (e  ),  which 

K  K  K  K 

indicates  a  reflected  wave  as  opposed  to  an  incident  wave  (see  Eqs.  2.40  to 
2.44)  and  that  e  must  be  replaced  by  e  .,  which  is  the  parameter  associated 

K  KJ 

with  the  reflected  wave. 
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APPENDIX  B 

THE  DEFINITION  OF  THE  FUNCTION  X  (e) 

IC 

The  complete  definition  of  the  multi-valued  function  ^(e)  is  of 
extreme  importance  in  the  solution  of  wave  propagation  problems  in 
orthotropic  media  by  the  method  of  self-similar  potentials.  The  solution 
of  the  algebraic  biquadratic  equation  (2.18)  yields  the  four  values 

(B.l) 

l3,e) . .  [bWwsn 

which  correspond  respectively  to  the  four  sheets  S-j ,  $2>  S3  and  of  a 
Riemann  surface  as  the  domain  of  definition  of  the  function  A^e)  [11]. 

The  four  sheets  are  connected  along  branch  cuts  in  a  manner  to  be  described 
later.  Since  the  ends  of  the  branch  cuts  are  the  singular  points  of  the 
function  A^Ce)  which  are  common  to  more  than  one  branch,  it  is  necessary 
to  determine  first  the  position  of  these  branch  points  on  the  Riemann 
surface.  Their  values  are  found  by  considering  the  following  system  of  two 
equations  [16]: 


If  x  f  0  then  x 

K  K 


p 

~+2bd9 "  and  back“substitution  t0  (B. 2a)  results  in 


Q(e)  =  0 


(B.4) 


The  polynomial  Q  (0)  (see  Eq.  2.20)  has  in  general  four  roots  which  assume 


real,  imaginary  or  complex  values  depending  on  whether  or  not  certain 
inequalities  among  the  elastic  parameters  are  satisfied.  The  results  of 


this  investigation  are  summarized  in  Table  B.l  where  it  can  be  seen  that 
there  exist  four  different  cases  denoted  by  I,  II,  III  and  IV. 

The  choice  of  the  branch  cuts  for  each  of  these  cases  is  shown  in 
Fig.  B.l.  Each  of  the  four  sheets  of  the  Riemann  surface  has  the  same 
branch  cuts  for  the  function  [(He^)]1*.  These  branch  cuts  connect  sheets 
S-j  to  S2  and  to  S^. 

The  remaining  branch  points  are  defined  by  the  condition  x^  =  0. 


Direct  substitution  into  Eq.  (B.2a)  yields  the  values 


The  location  of  these  points  on  the  Riemann  surface  depends  on  whether  the 
quantity  N2  is  positive  or  negative.  It  can  be  shown  that  when  N2  >  0  the 
branch  points  ©A  are  located  on  sheets  and  S3  while  ©D  are  on  sheets 
S2  and  S^.  The  choice  of  the  branch  cuts  in  this  case  is  shown  in 
Fig.  B.2  (a).  The  branch  cut  (-0^,  connects  sheet  to  S3  while 
(-0D»  0q)  connects  sheet  S2  to  S^.  When  N2  <  0,  both  pairs  +  ©A  and  +  0Q 
are  singular  points  common  to  and  sheets.  Figure  B.2  (b)  shows  the 
selection  of  the  branch  cuts  for  the  case  N2  <  0.  Again  the  branch  cuts 
(-0Q1,  -0D),(-©A,  ©A)  and  (©D>  0q^  )  connect  sheet  to  S3  while  the  branch 
cut  (“0qi>  0qi )  connects  sheet  S2  to  S^.  It  must  be  noted  here  that  N2  <  0 
is  not  compatible  with  the  inequalities  associated  with  cases  I  and  III 
(Table  B.l).  When  N2  <  0,  there  are  always  two  real  branch  points  + 
satisfying  the  condition  |0Q1|  _>  |l//d|. 

The  complete  Riemann  surface  is  given  in  Figs.  B.3  and  B.4  for  two 
distinct  cases.  The  construction  and  interconnection  of  the  four  sheets 
for  the  remaining  cases  is  made  in  a  similar  way.  Only  sheets  S-j  and  S2 
are  required  to  represent  the  mapping  of  the  wave  fields  as  discussed  in 
Chapter  4. 

The  important  behavior  of  *K(e)  for  values  of  0  on  the  real  axis  of 
the  Ime  _>  0  half  complex  plane  will  now  be  considered.  It  is  easy  to 


(B. 6d) 
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In  the  interval  l//a<  | Ree |  <  l//d,  A^(e)  is  imaginary  and  ( 0 )  is  real. 

For  any  values  with  | Ree |  _>  1//3"  the  behavior  depends  on  whether  N2  is 
positive  or  negative. 

When  N2  >  0 

\}{+  l//dj  =  (-N2/bd2)  55  (B.7a) 

which  is  imaginary  and 

a2(+  l//dT  =  0  (B. 7b) 

For  larger  values  of  Ree,  both  A-|(e)  and  X2 ( e )  remain  imaginary  unless  a 
real  branch  point  0O1  is  encountered,  in  which  case  they  become  complex. 

When  N2  <  0 

^(+1/40  =  0  (B.8a) 

and 

x2(±  l//d)  =  (-N2/bd2)  15  (B.8b) 

which  is  a  real  quantity.  The  existence  of  at  least  one  real  branch 
point  on  the  positive  real  axis  is  necessary  if  the  inequality  N2  <  0  is 
to  hold.  For  any  values  of  e  in  the  interval  l//cT<  |  Ree  |  £  Og-j ,  both 
A-j(e)  and  A2(e)  are  real.  For  larger  values  they  become  complex.  In 
Fig.  B.5  the  graphical  representation  of  the  real  values  of  A-|(e)  and 
X2 ( Q )  is  shown  for  real  values  of  e  in  the  indicated  intervals. 
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TABLE  B. 1 

THE  ROOTS  OF  THE  POLYNOMIAL  Q  (e) 


Conditions  of 

Elastic  Parameters 

Type  of  Roots 

Case  of  Cut 
Complex  Plane 

N, 

>0 

1 

Four  complex  (in  two 
conjugate  pairs) 

I 

N,  <  0 

O 

V 

CVJ 

Two  real  +  two 
imaginary 

II 

0 

A 

CVJ 

M  <  0 

Four  imaginary 

III 

M  >  0 

Four  real 

IV 

Note  1 .  The  roots  of  Q  (0)  are  given  by  the  values 


9  =  ±  Q01  =  ± 


9  -  ±  ®02  =  — 


M  -  /-4bdc\  y/2 


M  +  /  -4bdc2N]  ^1/2 

K1K2  J 


a)  Case  N2  >  0 


ImO 


b)  Case  N2  <  0 


FIG.  B.2  THE  BRANCH  POINTS  0.  AND  en,  AND  THE  CORRESPONDING  BRANCH  CUTS 


Section  T | 


Section  Tg 


Section  T3 


Section  T4 


o4 


FIG.  B.4  THE  COMPLETE  TYPE  II  RIEMANN  SURFACE  FOR  N2  <  0  (SECTIONS 

T-| ,  T2>  t3>  and  t4  indicate  the  interconnection  of  the 

FOUR  SHEETS) 
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APPENDIX  C 

THE  PARAMETRIC  EXAMINATION  OF  THE  GEOMETRY  OF  THE 
QUASI -TRANSVERSE  WAVE  FRONT 

The  quasi -transverse  wave  front  is  parametrized  by  real  values  of  in 
the  following  intervals 

a)  N^  >  0  -l//3"<_  l  l//d" 

b)  N£  <  0  ”®01  —  —  ®01 

1//T  <  |e]  |  _<  0Q1 

where  e2  represents  values  on  the  S2  sheet  and  e-j  values  on  the  S-|  sheet 
of  the  Riemann  surface.  The  two  cases  will  be  examined  separately  and 
only  for  the  positive  values  of  e^. 

C.l  The  Condition  N2  >  0 
It  can  be  seen  that 


(C.  la) 

(C.lb) 


4>2(0)  =  2b  >  0 

♦2(1//cT)  =  >  0 

while  from  Eqs.  (4.12) 

4>2(e2)  <0  for  N]  >  0 

4>2(e2)  >  0  for  N-j  <  0 


(C.  2) 


(C.3) 


As  a  result  1S  positive  throughout  the  interval  (0,  l//d)  and  thus, 
c»2  increases  monotonically  in  the  interval  (0,  v/2).  It  can  be  seen  also 
that 


tane2(02  =  0)  =  -  x2(0)  =  0 
tane2(e2  =  1//3)  =  -  x2(l//cf)  =  » 

i.e.  e2  varies  continuously  in  the  interval  (0,  u/2).  Since  x2 ( ©2 ) 
varies  monotonically  in  the  specified  interval,  the  behavior  of  B2  depends 
directly  on  the  function  il<2(02)  which  at  the  ends  of  the  interval 
(0,  l//d)  assumes  the  values 

2bN9 

*2<°>  =  •  W 

2  (C.4) 

<J'20  f/8)  =  -  2bdN^-  <0 


while  for  every  value  of  e2  in  the  same  interval 

8bdc2N, e9 

W‘73 


CQ°(e2)] 


(C.  5) 


It  can  be  seen  from  Eqs.  (C.4)  and  (C.5)  that  for  >  0  or  N-j  <0  and 
N3  >  0  the  function  ( e2 ^  remains  negative  throughout  the  interval  while 
for  N-|  <  0,  N3  <  0  changes  from  positive  to  negative  at 


fM-L  /? 


%  = 


N-|/ad 


K1 K2 


(C.  6) 


Using  the  equations  of  section  4.2  the  following  values  can  he  found. 
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4b2N, 

V°>  '  -  TT- 

02(1/^T)  =  -  l  020//d)]2  <  0 

0^(0)  =  0  (C*7) 

D^l/v'd)  =  0 
Let 

F(e)  =  (b-d)2  -  K^e4  (C.3) 

which  becomes  zero  at 

0F  =  (b-d)*5  /  (K1K2)!s  (C. 9) 

For  0  =  l//cf  the  function  F( e )  becomes  the  quantity  N4(see  Eq.  (2.1 0)). 

Let  first  N3  be  positive.  In  this  case  D2( ©2 )  remains  negative  in 
the  interval  (0,  1/i^d)  provided  one  of  the  following  three  cases  holds 

a)  N4  >  0 

b)  ^  <  0  and  N4  <  0  (C. 10) 

c)  N-j  >  0  and  D2(0p)  =  <  0 

As  a  result,  the  angle  B2  increases  monotonically  in  the  interval  (0.  ir/2) 
and  the  quasi-transverse  wave  front  is  a  strictly  convex  curve. 

If  conditions  >  0  and  >  0  are  satisfied,  then  D2 ( e2)  ^as  tw0 
roots  in  the  interval  (0,  l//cT)  denoted  by  0  and  0  .  The  values  of 
these  roots  can  be  conveniently  evaluated  by  using  a  numerical  solution. 
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1  *  9* 

It  is  not  difficult  to  see  then  that,  if  e  <  e  the  function  02(62) 

■j*  i*  2*% 

is  negative  in  the  interval  (0,  0  ),  positive  in  the  interval  (6,0) 

9* 

and  negative  in  the  interval  (0  ,  l/v^).  Therefore,  the  angle  62  increases 

l*  /  1*, 

monotomcally  from  0  to  6  (i.e.  the  value  of  ^  which  corresponds  to  0  ), 

1  ★  9* 

decreases  monotonically  from  e1  to  B  and  increases  monotonically  from 
2* 

6  to  tt/2.  The  quasi-transverse  wave  front  in  this  case  possesses  cusps 
in  each  quarter  plane. 

Consider  now  N3  <  0  which  implies  <  0.  In  such  a  case 
02(62)  changes  from  positive  to  negative  values  at  a  point  0  which  is  the 
root  of  D2(02)  in  the  specified  interval.  The  angle  B2  decreases 
monotonically  from  0  to  a  negative  value  -B  corresponding  to  0  and  then 
increases  monotonically  from  -6*  to  tt/2.  It  becomes  zero  for  a  value 
0  =0^  (Eq.  C.6)  which  is  the  root  of  function  ^2(02).  Thus,  the 

quasi-transverse  wave  front  has  cusps  located  adjacent  to  the  n-axis. 

C.2  The  Condition  N2  <  0 

The  functions  will  be  examined  in  the  interval  (0,  0^)  of  the  S2 
sheet  and  (©q-|,  l/>/d)  of  the  S1  sheet.  The  index  of  each  function  indicates 
the  appropriate  sheet  with  the  corresponding  interval.  It  can  be  seen  that 


♦,(l/«fl)  =  <  0 

*1 (eoi )  =  '  " 

♦2(0)  =  2b  >  0 


(C.ll) 


+2(e01)  =  + 


00 


while 


4>j  ( 6-| )  <0 

and 

*2  (02)  >  0 

Therefore,  «2  increases  monotonically  from  0  to  tt/2  for  a  continuous 
change  of  from  0  to  0q^  on  S2  sheet  and  from  0q^  to  l/*^3~on  S-j  sheet. 

Consider  the  following  quantities 
2bN, 

^2  ^  =  "  TwT 
h  (0oi) 

(C.12) 

*-,om  ■  -  >  o 

*  ■ 

while 

^2 ( 0 2 )  <  0  and  ^^ ( ©■] )  >0 

In  addition 

4b2N, 

D2<°>  =  -  SJ- 

D2(0O1)  =  '  “ 

(C.13) 

D1(0O1)  =  +  " 

D-j  (Vrf)  =  -  1  (l//?)]2  <  0 
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while 


0^(62 )<  0  and  Dj  (e-j )  >  0 

Let  >  0.  Then  functions  ^(e^  an£*  are  negative,  ^(e-i)  is 

positive  and  D^e^)  changes  from  positive  to  negative  at  a  root  e  .  The 

angle  32  increases  monotonically  from  0  to  tt/2,  continues  to  increase 

** 

monotonically  from  tt/2  toe  corresponding  to  the  value  e  and  then 
decreases  monotonically  from  that  value  to  tt/2.  The  quasi-transverse 
wave  front  has  cusps  adjacent  to  the  s-axis. 

When  <  0  the  behavior  of^e-i)  and  D-j ( B-j )  is  not  altered. 

However,  ^(e?)  and  D2(e2)  change  from  positive  to  negative  values  at 
the  root  e2*  =  (see  Eq.  (C . 6 ) )  for  ^2(e2)  and  at  e1  for  D2(e2).  As 
a  result  the  angle  32  decreases  monotonically  from  0  to  a  negative  value 
-3  (corresponding  to  e  ),  increases  monotonically  from  -3  to  3  (>  £) 

'2*  TT 

passing  through  32  =  0  at  e  and  62  =  -^  at  Qq-j  ,  and  then 
decreases  monotonically  from  3  (corresponding  to  e  )  to  a  value  of 
tt/2  at  e-  =  lA'cT.  The  quasi-transverse  wave  front  possesses  cusps 
adjacent  to  both  5  and  n  axis.  A  summary  of  this  investigation  is 
given  in  Table  1 . 


appendix  d 


NUMERICAL  TECHNIQUES 


D.l  The  Square  Root  of  Q  (0  ) 

The  polynomial  Q  (e^)  given  in  Eq.  (2.20)  can  be  written  in  the 


q  <<g  -  k,k2  (eK-0ol)  (eK+e01)  <VV  <V<te> 


where  +  and  ±0O2  are  its  four  roots.  Let 


V8oi  ■  rn  e 


Ve0l  =  rI2  e 


V902  '  r21  e 


V0O2  =  r22  e 


The  square  root  of  Q(e)  is  then  qiven  by  the  expression 

q  <<g  =  [Q(«k)11/2  ■  ±  VT (9k-O0,)1/2  (VV1 
x  <  <VS02 >1/2  ■  iVT  rq  e'*q 


where 


rq  =  trirr12'r21'r22] 


and 
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_  <f’ll  +  *12+^21  +  ^22  ,r 

The  range  of  the  arguments  (j^,  <j>22>  depends  on  the  relative 

position  of  the  roots  +0g-|  and  i^and  is  different  for  each  of  the 
four  types  of  Riemann  surface.  The  sign  in  Eq.  (D.3)  is  fixed  for  each 
case  in  such  a  manner  as  to  yield  a  positive  value  of  q(e  )  for  e  =  0. 
Figure  D.l  shows  the  polar  rays  and  arguments  to  determine  q ( e ^ )  on  a 
complex  plane  with  the  branch  cuts  of  a  Type  II  Riemann  surface.  Note 
that  in  this  particular  case  K2  <  0,  therefore,  =  i  /j K-j  K2j.  In 

accordance  with  the  previously  mentioned  condition,  Eq.  (D.3)  must  be 
chosen  with  the  minus  sign. 

D.2  The  Square  Root  of  A  (e  ) 

_ _ _ K  K 

For  every  value  of  the  variable  e^,  the  quantity 


A  (e  ) 

K  K 


b+d-lV  +  ( -  T )  Kq  ( e  ^ ) 
2bd 


(D.5) 


assumes  values  in  a  complex  domain  S  .  The  function  q(e  )  is  given  in 

Ak  k 

Appendix  D.l.  The  complex  domain  S ^  has  a  branch  cut  along  the  real 
A  axis  from  point  0  to  infinity.  It  is  convenient  to  choose  this 

K 

branch  cut  in  such  a  way  that  for  all  values  of  in  the  upper  half 
S  sheet  (i.e.  Ime  >  0)  the  values  of  A  (e  )  lie  on  a  single  sheet  of 

K  K  <  < 

the  S  domain.  It  turns  out  that  for  certain  cases  the  branch  cut  may 
A 

K 

most  conveniently  be  taken  along  the  positive  real  axis  while  in  the 


remaining  cases  it  is  to  be  taken  along  the  negative  real  ax7S  of  the 
domain  .  Once  this  choice  has  been  made  appropriately,  the  evaluation 

K 

of  the  function 


A  (e  )  =  A  (e  ) 

K  K  K  K 


is  a  simple  matter.  The  sign  in  Eq.  (D.6)  has  been  fixed  from  the 
mapping  procedure  which  is  presented  in  Chapter  4. 


D.3  The  Evaluation  of  e  near  the  Wave  Fronts 


In  section  5.2.2  it  has  been  shown  that  for  points  near  the  wave 
fronts  the  rate  of  change  of  with  time  is  very  large.  The  increased 
curvature  of  the  surface  that  is  formed  from  the  values  of  the  right  hand 
side  of  Eq.  (5.1)  in  this  range  of  values  of  calls  for  a  more  accurate 
evaluation  of  the  predicted  value  in  the  predictor-corrector  iterative 
method  than  the  one  obtained  from  the  utilization  of  the  trapezoidal  rule. 
This  is  achieved  by  using  an  integration  formula  which  is  obtained  as 


follows. 


From  the  discussion  of  section  5.2.2  it  can  be  written  that 


V9k0  5 


Using  this  equation  as  well  as  Eq.  (5.3)  the  differential  equation 
given  in  Eq.  (5.1)  can  be  expressed  approximately  as 


[A  (e  )+(t-tn)  B(e  )] 


ICO 


The  result  of  the  integration  of  Eq.  ( D. 8)  between  tffl  =  tp  +  m  At  and  t^ 
t0+(m+l)At  is  written  by  the  integration  formula  [33,  Appendix  A] 

!] 


A6 


=  SEi  |gm  j^2(l+m)  [(m+l)1/2  -  m1/2]-  |  [(nH-1 ) 

*  v,  [!M>3/2  -  -  J'2]' 


3/2_m3/2 | 


Here,  m  is  an  integer  and 


g  =  / t  -tn  A ' (e  )  =  A(e  )+B(e  )  (t  -tn) 
3m  m  0  k  '  k  2  k  k  m  0 


(D.9) 


Vi  =  A<eK>+BteK>  (Vr‘o> 


(D.10) 


where  e  is  the  known  value  at  time  t  . 

<  m 

The  procedure  to  be  followed  in  applying  the  predictor-corrector 
method  in  its  modified  form  is  now  briefly  explained.  Let  the  known  value 
of  0^  at  time  t^  be  written  as  A  first  approximation  at  time  tffl+1 

is  given  by 


e1  =  e  +  Ae1 
k,  m+i  K,m  k,  m 


(D. 11 ) 


where  Ae1^  m  is  obtained  from  Eq.  (D.9)  and  g^,  g^+1  are  given  by  Eq. 
(D.10).  Subsequent  iterations  yield 


e  =  e  +  A8 

K,m+!  K,m  K,m 


(D.12) 


where  Ae  is  evaluated  from  the  integration  formula  given  in  Eq.  (D.9) 

K  y  Ml 

and 


g_  =  /t-tn  (e  )  =  gm 
m  m  0  k  K,m'  3m 


Vl  =  /  tm+l"t0  6k  (0K’1m+l) 


(D.13) 
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At  m=0  the  procedure  requires  special  attention.  It  can  be  shown  that 


so  =  A(e«<P  * 


K<6*0> 


A  first  approximation  of  the  value  of  at  t^  can  be  obtained  by 
integrating  Eq.  (5.1)  written  as 


89« 

31  ^ 


(D. 14) 


(D. 1 5) 


between  tg  and  t-j ,  i.e. 


e1  .  =  e  _  +  ag1  n 

K,1  <0  <,0 


(D. 16) 


where 


0  =  2  /tl't0  A 


(D. 17) 


.2 


The  second  of  Eq.  (D.13)  can  now  be  used  to  obtain  g^  which  when 
substituted  in  Eq.  (D.9)  together  with  the  value  from  Eq.  (D.14)  yields  a 
new  value  of  Ae  .and  so  on. 

K 

The  advantage  of  this  numerical  technique  lies  on  the  fact  that  it 
produces  a  very  good  first  approximation  resulting  in  very  fast 
convergence. 

D.4  Integration  Formulas  for  Points  Near  the  Ends  of  the  Contours 


Three  integration  formulas  will  be  given.  Which  of  these  should  be 
applied  depends  on  whether  the  point  under  consideration  lies  on  a  wave 
front  or  not  and  on  the  nature  of  the  singularity  caused  by  the  expression 
G< ( e< ) •  ™ese  are  cases  in  which  a  standard  numerical  integration  technique, 
such  as  Simpson's  rule,  cannot  be  directly  applied. 


1   mmmmami 


D.4.1  First  Integration  Formula 

Consider  the  integrands  of  Eq.  (3.36)  which  contain  the  function 
6  (e  )  in  the  denominator.  Assume  also  that  the  end  values  of  the  contour 

K  K 

C  do  not  represent  points  on  the  wave  fronts.  Then,  these  integrands  at 

K 

values  of  e  very  close  to  the  ends  of  the  contour  C  can  be  approximated  by 

K  K 

expressions  of  the  form 

A„+B„(8  -e")  F(e")  . 


F  (eK)  - 


-  n  n  <  k  1  . 


(D. 18) 


/e  "  0  n 

K  tcO 


V6ko 


where  0kQ  represents  the  end  of  the  contour  and  the  beginning  of  the 


integration  interval.  The  constants  An  and  Bn  are  to  be  determined.  The 
contribution  of  the  integral  between  values  o'1  and  en+^  can  now  be  written 

K  K 
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n 


0.4. 2  Second  Integration  Formula 

This  formula  is  used  with  the  integrand  expressions  that  contain  the 


function  G^(0^)  in  the  numerator.  It  is  applicable  for  points  not  on  the 


wave  fronts.  In  this  case  the  integrand  is  approximated  by  the  expression 


F(eJ  =  /0."e.n  =  F(eJ  /07"0 


K  < 0 


V  <0 


(0.21) 


and  the  contribution  of  the  integral  between  and  is  given  by 


CWe^-el)]  /T^  deK  = 


n  n  <  tc 


2B 


*  -Vs  Cte^’-e  0S/2  -  (e"-e  n>5«] 


<0 


(0.22) 


Here,  the  constants  A  and  B  can  be  determined  by  Eq.  (0.20)  for  each 

n  n 


interval  of  integration. 


D.4.3  Third  Integration  Formula 

This  is  a  variation  of  the  first  integration  formula  for  points  on 
the  wave  fronts.  In  such  a  case  the  function  ^(e^)  is  expressed  as  in 
Eq.  (5.9  ).  The  same  integrands  that  were  given  by  Eq.  (D.13)  can  now 
be  written  as 


F(e  )  =  WV6*)  __ 


(0.23) 


VeK0 


0  -0  n 
K  KU 


The  result  of  the  integration  between  0n  and  0n+^  is 
0n+l 


J3  = 


<  VBn<VO 


n 

,0K 


0  -0  n  de* 

K  KO 


-■i..  ...  V 


.*■  -k-"  .1. 


FIG.  D.l  THE  ARGUMENTS  FOR  THE  DEFINITION  OF  q=[Q(e  )]1/2 
ON  THE  TYPE  II  RIEMANN  SURFACE  K 
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APPENDIX  E 

THE  VELOCITIES  OF  THE  WAVES 


In  anisotropic  media  the  normal  nK  to  the  tangent  plane  at  a  point 
P  of  the  wave  front  is  not  parallel  to  the  ray  at  the  same  point 
(Fig.  2).  The  normal  n  defines  the  velocity  vector  while  the  ray 
defines  the  path  of  energy  flux.  In  orthotropic  materials  the  two 
directions  coincide  only  along  the  axes  of  elastic  symmetry. 

Let  l  and  m  be  the  direction  cosines  of  the  vector  n^.  The 
equation  which  determines  the  wave  velocities  associated  with  n^ 
is  given  from  the  determinant  [21,  5208] 


2  2  2 
a£  +dm  -u  c 
< 

dm 


dm 

de2+bm2-u  2 
< 


=  0 


which  results  in  the  biquadratic  equation 


ok  - r(a+d)£2+(b+d)m2]  v2+[al2+dm2)  (d£2+bm2)-c2£2m2=0 


(E.l) 


(E.  2) 


The  solution  can  be  written  as 

2 _ ] _ 

K  e2  +  x2(e  ) 

K  K  K 


(<=i,2) 


(E.3) 


Every  point  P  of  the  wave  fronts  is  uniquely  determined  by  a  real  value  of 
the  parameter  e  .  By  Eq.  (E.3)  this  value  determines  the  magnitude  of 
the  velocity  of  propagation  of  the  corresponding  wave  in  the  direction 
specified  by  nK. 

In  transversely  isotropic  media  the  velocities  associated  with  the 


187 


vector  R  which  has  direction  cosines 
determinant 


2 

u 


YU  -  o 

y12 

2 

y12 

y22  '  u 

y13 

y23 

where 

2  2? 
Y-j  i  =  a  t  +dm  +en 

y22  =  bm2+d(£2+n2) 

2  2  2 
Y33  =  +dm  +an 

>12  ’  <*» 
y13  =  (a-e)  In 

y23  =  cmn 


m  and  n  are  given  by  the 

Y13 

y23  =  0  (E.4) 

2 

Y33  “  u 

(E.5) 


When  the  vector  N  lies  in  the  (x,  y)  plane  (i.e.  n  =  0),  Eq.  (E.5)  is 
equivalent  to  Eq.  (E.3)  which  gives  the  velocities  of  the  quasi¬ 
longitudinal  and  quasi-transverse  wave,  and  to  the  equation 

u2  =  el2+dmZ  (E.6) 

which  determines  the  velocity  of  the  third  wave. 
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APPENDIX  F 

SPECIAL  RELATIONS 

In  the  development  of  the  solutions  of  the  present  study,  the 
algebraic  work  may  be  reduced  by  the  observation  that  certain  functions 
can  be  written  in  equivalent  forms  which  are  more  convenient  to 


manipulate.  Some  of  these  relations  are  given  here; 

x|(  e)+x|(  e)  =  (b+d-Le2)/bd  (F-1) 

x^(e)-x|(e)  =  -4TeT/bd  (F.2) 

A2(e)  A2(e)  =  (1-ae2)  (l-de2)/bd  (F.3) 

[ae2  +  dx*(e)-l]  [ae2+dx2_K(e)-l ]  =  c2e2(ae2-l)/b  (F.4) 

[ae2  +  dx2(e)-l]  [de2  +  bx2(e)-l]  =  c2e2  A2(e)  (F.5) 

[dx^  ( e)-l ]  [ae2+dA2_K(e)-l]  =  e2[ad(de2+bA2(e)-l}-c2]/b  (F.6) 

d[ae2+bA2(e)-l]  [de2+bA2(e)-l ]  = 

=  (b-d)  (1-ae2)  [de2+bA2(e)-l ]  +  bc2e2A2(e)  (F.7) 

K  K 

K1K2  =  l2  ‘  4ab<j2  (F.8) 

L  (b+d)  -  M  =  2bd  (a+d)  (F.9) 

M2+4bdc2N1  =  (b-d)2  (F.10) 

L  (b-d)  +  M  =  2bN3  (F.ll) 

L  (b-d)  -  M  =  2d  [d(b-d)  +  c2]  (F.12) 


Is  given  by  the  expression 
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APPENDIX  G 

ISOTROPY  AS  A  SPECIAL  CASE 

An  anisotropic  material  of  the  class  considered  in  this  study 
which  has  elastic  parameters  satisfying  the  conditions 

a=b 


c*a-d 


(G.l) 


e=d 


forms  the  special  case  of  isotropy.  Let  X  and  y  be  the  Lame  constants. 
Then,  it  is  net  difficult  to  see  that 


X  =  p(a-2d) 
y  =  pd 


(G.2) 


Whenever  Eqs.  (G.l)  hold,  the  expressions  of  Chapters  2  and  3  yield  the 
solutions  to  dynamic  problems  for  an  isotropic  material  with  constants 
defined  by  Eq.  (G.2).  It  is  useful  to  give  some  of  the  expressions 
which  were  used  to  develop  the  anisotropic  solutions  in  their  simplified 
form  for  isotropy,  i.e. 


K1  =  4d  (a-d) 

K2  =  M  =  N1  =  0 
N2  =  N3  =  d(a-d) 
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L  =  2ad 

Q  (e)  -  (a-d)2 

x-j  ( e )  =  (a"1  -  e2)1* 

A2(e)  =  A3(e)  =  (d_1  -  e2)*5 

^(e)-»|(e)  -  - 

ae2+dA2(e)-l  =  -  (a-d) (a_1-e2) 

de2+bA2(e)-l  =  -  (a-d)e2 
ae2+dA2(e)  -1  =  (a-d)e2 

de2+bA2(e)  -1  =  (a-d)  (d-1  -  e2) 

[x1(e)-x2(e)3  R(e)  =  d(a-d)  (a"1  -  e2)5  Rj(e) 

where 

Rj(e)  =  (d"1  -  2e2)  +  4e2(a-1  -  e2)*2  (d-1  -  e2)*2 
is  the  Rayleigh  function  for  isotropic  media. 
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